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THE SKYLAB CONCENTRATED ATMOSPHERIC RADIATION PROJECT 
CONTRACT T-4714B 


P. M. Kuhn, W. E. Marlatt, and V. S. Whitehead 


In conjunction with SKYLAB measurements, data were obtained at 
Earth's surface, in the boundary layer, and from aircraft, to be 
used to improve understanding of radiation transfer within the 
atmosphere. Measurements included total net radiation, albedo 
and temperature of the interface, atmospheric profiling of atmos- 
pheric radiative flux and layer cooling, and aerosol density, 
size, and distribution; total extinction coefficients were 
calculated. Instrumentation for observations is described. 

Data from radiometersonde observations, and the computer listings 
of the data are included in Appendices. Models for computing 
atmospheric IR transmission are compared. Spectral estimates of 
planet albedo are given and compared with SKYLAB observations. 
Performance of models is evaluated, and recommendations are made 
regarding the need for further studies. 


1. INTRODUCTION 

Participants in the SKYLAB Concentrated Atmospheric Radiation Project 
(SCARP) utilizing the Earth Resources Experiment Package flown on SKYLAB 
were Dr. P. M. Kuhn of the NOAA, Environmental Research Laboratories, 

Dr. W. E. Marlatt of the Department of Earth Resources, Colorado State 
University, and Dr. V. S. Whitehead of NASA, Earth Observations Division of 
Johnson Space Center. These were joined in the Phoenix field measurements 
by R. F. Pueschel of the NOAA Environmental Research Laboratories and 
L. Edwin Williamson of the White Sands Missile Range, New Mexico. 


1.1 Approach 

The ultimate purpose of the project was to arrive at a more complete 
understanding of radiation transfer within the atmosphere, including the 
contribution of aerosols to this transfer. Intermediate goals were as 
follows: 

1. To acquire data giving a comprehensive description of the 
atmosphere's structure and composition over sites of differing 
radiation properties for a variety of air masses. Our measure- 
ments were taken at the surface and at levels between the surface 
and the spacecraft to support the SKYLAB infrared spectrometer, 
and scanner measurements of selected sites. 

2. To acquire and implement several models of radiation transfer 
for both visible and thermal regions; and to test these models 
against observations. 






3. To apply the result of these tests to improve models for the 
variety of existing air masses and atmospheres. 

Data were acquired in two related experiments, in which a variety 
of measurements were taken (1) at the Earth's surface and in the boundary 
layer, and (2) from aircraft. A number of existing radiative trahsfer 
models allowed for varying air masses and amounts of wet and dry aerosols. 

It was anticipated that other investigations could make use of the data 
base acquired. Hence the experiment was declared an open one and other 
investigators participated on an exchange-of-data basis. Among them were 
Ed Williamson of White Sands Missile Range and Rudolf Pueschel of NOAA's 
Atmospheric Physics and Chemistry Laboratory. Data exchanges occurred with 
EPA which provided a lidar in the first field effort, with the Texas Air 
Control Board for data acquired in the Houston area, and with other SKYLAB 
investigators. 

To acquire the desired variety in the data base, sites were selected and 
data acquired onshore and offshore in the vicinity of Houston, Texas, to 
provide a warm moist clean and a warm moist dirty atmosphere; over White 
Sands Missile Range, N. M. , (both sarfds and the mal pais) to provide extremes 
of surface characteristics for a clean dry atmosphere; and over Phoenix, 
Arizona, to provide hot dry clean and hot dry dirty atmosphere. Other sites 
surveyed (San Francisco, Four-Corners, and St. Louis) were not used because 
of operational constraints, primarily those of weather. 


1.2 Acquisition of Data 

Each field effort differed somewhat from the other in the type of data 
acquired as the investigators attempted to apply the experience of the 
preceding effort. Description of a typical exercise follows. 

1. Personnel and equipment were moved to facilities nearest to the 
points of the field effort about three days prior to the scheduled 
SKYLAB overflight. Last minute changes to the operation plans were 
made. 

2. Two days prior to the SKYLAB overflight a low key field effort was 
initiated primarily to check out equipment and procedures. During 
this day one or more radiometersonde releases would be made by 

P. M. Kuhn and Associates. Surface meteorological and radiometric 
measurements would be made at more than one site. The Colorado 
State aircraft would fly at least one multiple level pattern 
between sites with sensor operating for general familiarity, timing, 
and instrument test. This days' effort would be concluded by a 
general debriefing on problems and re-planning where desirable. 

3. On the day prior to the SKYLAB overpass, a schedule was followed 
identical to that planned for the day of the overpass. The 
schedule included hourly, or continuous measurement of the 
following: 
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(1) Downwelling and upwelling solar flux. 

(2) Direct beam solar radiation and albedo. 

(3) Total net radiation at the surface. 

(4) Radiation temperature of the surface, sky, and cloud. 

(5) Wind, temperaturt:, and dew point. 

(6) Sky cover (photographs). 

(7) Radiometersonde profiles (temperature, dew point, up- 
and downwelling thermal radiation, and layer cooling) 
each three hours. 

(8) Interface temperature measurements from ground and 
helicopter. 

(9) Aircraft profiles from near surface to above the haze 

layer (profiles of aerosol number density and size 
distribution, upwelling and downwelling solar flux, 
and. profiles of temperature and dew point). Three 
times daily: early morning, at time of day of SKYLAB 

overpass and late afternoon. If the time of day of 
SKYLAB overpass differed greatly from the time of 
maximum surface heating another flight would be 
attempted at that time, 

4. On the day of the overpass, the same schedule was followed, with 
the addition of the SKYLAB crew's participation. Because the 
SKYLAB data had to be acquired within 30 seconds, the crew had 
been provided photographs of measurement sites, to aid in recog- 
nition and fast response when the spacecraft passed over the site. 
SKYLAB data were obtained for us with an S190A multi spectral 
camera, an S190B mapping camera, an SI 91 spectrometer, and an 
S192 mul tispectral scanner. The camera and scanner system were 
pointed at nadir in the Z-local vertical mode. The spectrometer, 
however, could be directed to point at a specific site by the 
operator to acquire several spectral scans. It should be noted 
that late in the planning stages for SCARP, the opportunity to use 
a helicopter mounted SI 91 spectrometer was accepted and profiles 
from near-surface to 12,000 feet were acquired for the same sites 
as those observed by SKYLAB. 

5. It was initially planned to maintain the same observation schedule 
the day following the SKYLAB overpass as on the two days preceding 
in order to take advantage of the established field deployment to 
acquire data from more varied air masses. This was discontinued 
about midway through the field program for efficiency in scheduling 
Exercises were then terminated about 24 hours after the SKYLAB 
overpass with only a single, or at most, two flights by the 
Colorado State University aircraft and one or two radiometersonde 
releases on that day. The exercise was concluded with a review 

of the effort and suggestions for planning of the next field 
exercise. 
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2. . SURFACE AND BOUNDARY LAYER EXPERIMENT (NOAA) 
P. M. Kuhn and L. P. Stearns 


2.1 Instrumentation 


2.1.1 Precision Radiometer 


The Barnes PRT-5 radiometer was used for measurements of the interface 
temperature at the various SKYLAB sites. It consists of an optical unit and 
an electronic unit with interconnecting cables. The optical unit compares 
the amount of energy emitted by the interface with that emitted by an inter- 
nal, controlled reference. The electronics unit converts the comparison into 
a voltage that can be recorded or read on a front panel meter. The detector 
is a hyper- immersed thermistor bolometer in a reference temperature cavity. 
Lenses restrict the spectral interval measurements from 8-14 microns. The 
instrument accuracy is 0.5C and the response time varies with the setting of 
an adjustable bandwidth setting from 5 to 500 milliseconds. The field ol 
view is approximately 2“. 


Effective radiance (N^^^) is defined as 

rv, R.. 


^eff = / ^ 
0 ^pk 


NgB dx, 


( 1 ) 


R. 


when p — is the normalized spectral response of the instrument and Ngg is the 
black body radiance. 


2.1.2 The Radiometersonde 


The radiometersonde used for soundings of the atmospheric temperature, 
pressure, relative humidity, and hemispheric radiation flux profiles, looking 
upward and downward at the various SKYLAB sites, was a combination of the 
Suomi -Kuhn Economical Net Radiometer attached "piggy back" to a standard 
balloon-borne radiosonde. The radiometer is a double-faced hemispheric 
bolometer with broad-response blackened' sensing surfaces at night, white 
in the sunlight, shielded by two thin polyethylene membranes. The sensors 
used were rod thermistors with precalibrated constants. Flux evaluations 
are made for 0.5-micron spectral intervals to compensate for the irregularity 
of the blackened surfaces. The unit is attached to a radiosonde in a 
horizontal manner with a sequencing device enabling the data from the radi- 
ometer to be recorded in sequence with that of the radiosonde. These data 
from the balloon-carried radiometersonde are transmitted by a 403-MHz trans- 
mitter to a ground receiving location and recorded on a strip chart. At the 
White Sands and Phoenix sites, the transmitter was a 1680 MHz unit. 


The frequencies representing temperatures were recorded and changed to 
temperatures based on the expression 


( 2 ) 


F X R £ 
f = ref ref 

^therm ^ *^ref 

where is the resistance of the modulator equal to 44.9 kil ohms and 

is the frequency of the modulator at this resistance equal to 190 cps 
or 95.0 frequencies on the strip chart record. 

The flux at any pressure level, looking upward or downward, may be 
obtained by the equations 

F4- = aTj + A[-k.(Tjj-T^)/D - k^(TQ-T^)/d] + Ad T^/dt, 

Ft = oTjj + A[+k.(Tjj-T^)/D - k,^{Ta-Tj^)/d] + Ad T^/dt, 

and = Ft - F+. 

Since A, B, k^ , k^, kj^, a. A, D, and d are all defined constants, it 
remains only to measure the parameters T^, (top plate temperature), and 
T|^ (bottom plate temperature) to determine the radiant fluxes with the 

predetermined log term based on the time, t. Once the net radiative flux 
is calculated, an additional term may be computed, that of the uncompensated 
atmospheric layer warming/cooling (W). It may be expressed as 

g (Fi-F, i) 

w = - IMO , (3) 

where g is the acceleration of gravity, Cp is the specific heat at constant 

pressure, and 1440 is the number of minutes in a day. The error of the 

instrument is 7 watts per square meter. 

2.1.3 The Net Radiometer 

The Suomi -Franssila-Islitzer Net Radiometer was used for measuring the 
net flux of radiation energy through a surface parallel to the Earth's surface 
at the SKYLAB sites. This instrument has a specially designed vane in the 
nozzle throat and an electric heater to furnish a sensitive control over the 
ventilation so that the cooling power in each side of the plate can be 
equalized. The sensing range is from 0.17 to 80 microns from a sensor con- 
sisting of a glass plate wound by 360 spaced turns of copper and constantan 
wire to form a thermopile after proper plating and blackening. About six 
percent of the entire surface of the blackened sensor is also coated with 
a strip of highly reflective white paint to compensate for solar radiation. 

The resistance of the thermopile is read directly from the calibrated sensor 
and is recorded in millivolts. It is capable of measuring total net radiation 
to within an accuracy of 2 percent. 
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2.1.4 The Pyronometer 

The Model 2 Eppley Pyronometer was used to measure the incoming solar 
radiation and albedo at the SKYLAB ground sites. This is a precision 
pyronometer incorporating a fast-response, wi rewound, plated with copper/ 
constantan, temperature-compensated thermopile and two miniaturized hemi- 
spheres of W6 7 glass. It was calibrated against standards verified against 
the Eppley primary reference group (Eppley-Angstrbm electrical compensation 
pyrheliometers) which maintains and reproduces the International Pyrheli- 
ometric Scale in the United States. 


2.2 Radiometersonde Observations , 

During SKYLAB Phases II, III, and IV, 28 radiometersonde ascents were 
made before, after, and during the overpass times of the satellite laboratory 
They occurred in the oil fields near Rosenberg, Texas; at White Sands, New 
Mexico, at both the lava area and the sands area; in the desert near Phoenix 
Arizona; and finally in St. Louis, Missouri. At the last site, no SKYLAB 
measurements were made, but our data are included. Table 1 is a list of the 
ascents with comments. 

2.2.1 The Data 

In Appendix A are the data plotted by a CRT unit of CDC 6600 of the 
ascents. Air temperature is denoted by T, upward radiant flux by A, down- 
ward radiative flux by v, net radiation flux by 0, warming/cooling by C, 
and humidity by M. Here langleys/min may be rewritten cal cm'^min"^. The 
computer plots the humidity in gm kg-i, using the radiation scale. 

Appendix B is the computer listing of the data in Appendix A. Column 
headings are provided and again ly/min is equivalent to cal cm'^min-^. 

The relative humidity is in percent. In flight number one, only, the 
radiative flux profiles are invalid and only the pressure, temperature, and 
moisture profiles should be considered. 


2.3 Thermal Radiation Calculation Model, Program Radiance 

The infrared radiation transfer program RADIANCE has been developed and 
modified to include ten options, all involving use of the radiative transfer 
equation. The program takes into account transfer through aerosols, the 
latest measurements of carbon dioxide, water vapor, window regions, and 
ozone. Therefore, radiance and irradiance may be determined for any infrared 
spectral interval to one inverse centimeter and at any angle either to 
include or to exclude aerosols or carbon dioxide. It may calculate black- 
body equilibrium temperatures, spectral intervals, weighting functions, 
brightness (apparent) temperatures, net irradiance, and profiles, or it may 
be inverted to infer water vapor .at specified levels when only the tempera- 
ture profile is known. Appendix C is a logic flow chart of the program and 
Appendix D is a listing of the program. 
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Table 1. Badiome ter sonde Ascents 


No. 

Date 

Location 

Time 

Weather/Clouds 

1 

6/04/73 

RNB 

1300 CDT 


3 

6/05/73 


0537 CDT 


5 

6/05/73 


1713 CDT 


6 

6/06/73 


0855 CDT 

5Cu 

7 

8/07/73 

RNB 

1715 UT 

5Cu 2Ac 2 Ci 

8 

8/07/73 


2212 UT 

35 <f) Cu 7 Ci 

9 

8/08/73 


0146 UT 

9 Ci 

10 

8/08/73 


1438 UT 

2 Cu 30 , 1 Ac 120, 1 Ci 

1 1 

8/08/73 


1600 UT 

4 Cu 

12 

8/08/73 


1937 UT 

4 Cu 30 , 4 Ac 120, 2 Ci 

13 

8/09/73 


1516 UT 


li» 

8/09/73 

BUC TWR 

0835 CDT 


15 

8/11/73 

W.S. LAVA 

2138 UT 


16 

8/11/73 


0000 UT 


17 

8/12/73 


1432 UT 


18 

8/11/73 

W.S. DUNES 

2130 UT 


19 

8/11/73 


2400 UT 


20 

8/12/73 


1437 UT 


21 

9/06/73 

PHX #5 

0800 MST 

0 Ac, Ci 

22 

9/06/73 


1015 MST 

0 Ac, Ci 

23 

9/06/73 


1405 MST 

0 Ac, Ci 

24 

9/06/73 


1900 MST 

0 Ac, Ci 

25 

9/07/73 


0800 MST 

Clear 

26 

9/07/73 


11 30 MST 

Clear 

27 

1/18/74 

ST. LOUIS 

1250 CDT 

5 Cc, Cs 

28 

1/18/74 


1535 CDT 

E5, 10 Cu RW - 


Comments 


Air temps only 


N i te 
Overpass 


Thru front 


7 




1 


2.3.1 Purposes of RADIANCE Program Sections 

Radiance and irradiance . To determine the radiance and integrated 
irradiance values over desired wavelength intervals, angles and specified 
filters with input pressures, temperatures, and mixing ratios. This also 
provides the integral water vapor and the carbon dioxide and water vapor 
totals above or below the background pressure. 

r 

Black- body equilibrium temperatures . To determine the black-body 
radiance and normalized radiance for specified frequencies and filters over 
an input temperature range at specif'jed intervals. 

Spectral intervals . To list the radiance of each frequency of a band 
at ten inverse centimeter intervals. 

Infrared weighting . To determine the quotient of the differential of 
the transmissivity for two levels divided by the differential of the natural 
logarithm of the pressure levels printed for each average pressure and 
temperature in a temperature versus pressure profile for each specified 
frequency and angle. Quotients are radiometric transmission weighting 
functions enabling one to determine "penetration" depth for a radiation 
observation in the atmosphere. Output includes punched cards of average 
pressures and quotients. 

Atmospheric corrections to observed brightness temperatures . To adjust 
the apparent or brightness temperature of the background or interface, for 
specified intervals and angles. This procedure reconstructs the physical 
temperature of the interface. 

Biqnell method of determining aerosol contribution . To determine 
radiance and integrated irradiance using two absorption coefficients in the 
continuum at specific wavelengths following the method of Bignell (1970). 

Aerosol . To determine radiance and integrated irradiance for specified 
layers in an atmospheric profile for different types and sizes of aerosols, 
using a bulk volume absorption coefficient. 

Net integrated irradiance and layer warminq/coolinq . To determine the 
difference between the upward and downward irradiance, using the differ- 
ential irradiance to determine the warming/cooling of a layer. Output 
includes punched cards of pressure, temperature, and layer warming/cooling. 

Water vapor inference . To determine the total water vapor burden above 
the reference level and the mixing ratio at the reference level by inference 
of the radiative transfer equation. To determine the relative humidity at 
the reference level with respect to water. 

Profiler . To determine a temperature or humidity profile from aircraft 
by measuring the radiance with assorted carefully chosen filters or at 
different angles by inverting the radiative transfer equation and using 
the already determined infrared weighting functions. Output includes cathode 
ray tube plot. 
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2.3.2 Method of Development of the Transfer Model 

Determination of the radiance . The transfer of radiation observed or 
calculated at a reference level, r, in the atmosphere is given by 

= - J J (j){v) B(v,Tg) dx dv + f B(v,T )(|) (v)x dv, (4) 

Vi Vi 

where v is the wavenumber, 

T is the atmospheric transmission in specified spectral intervals 
and, by the multiplicity theorem, is equivalent to 

’^H20* ^CO^* ^03 * ^conti nuum, 

(f) is the response function of the radiometer including the filter 
and detector transmission product, 

B is the Planck function, 

T^ is absolute temperature of an air layer, 

Tq is the physical temperature at the surface, and 

dx is equivalent to x.^^ - x., where i refers to successive 
atmospheric levels.' ' 

With no atmosphere present Eq. (1) may be written. 


^2 

No = / B(v,Tq)c))(v) dv. (5) 

The Planck body radiation expression is given by 

■^2 

j B(v,T) dv = I (av3/exp(bv/T)-l)Av (6) 

Vi V 

where a is equal to 3.7413x10-5 erg cm2sec-i and b is equal to 1.4389 cm deg. 

Combining Eqs. (2) and (3) permits the additional solution of T, the 
equivalent black-body or brightness (apparent temperature) given by 


T = I bv/ln ((^(v)vVNo)+1)Av. (7) 

V 

Where, in the usual sense, an atmosphere is present, Nq must be replaced by 
Nr from Eq. (1) and the solution proceeds to T, a brightness or apparent 
temperature not equal to the physical temperature Tq given in Eq. ^1). "T" 

is the required temperature of any interface. The brightness temperature 
T, is adjusted to give the true, physical temperature, To, by radiometric 
vertical profiling of the interface temperature through various atmospheres. 
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Addition of aerosols to model , (a) Bignell continuum transfer model : 
Water vapor continuum is determi ned as suggested by Bignell (1970) by using 
the total absorption coefficient at temperature T and total pressure p (p,e) 
by 

k(T,p,e) = l<i(t) p+k 2 (t)e (8) 

where ki is the absorption coefficient at unit total pressure for foreign 
broadening and kz is the absorption coefficient at uhit partial vapor 
pressure for e-type absorption, (b) Particulates transfer modeling: The 
transmissivity of aerosols may be determined by 

-In T = K / dz (9) 

where K is the volume absorption coefficient, and z is the height of the 
layer specified. 

Layer warming/cooling . It has been determined by Robinson (1950) that 
the use of radiance measurements at 52.5 degrees off zenith in calculating 
the irradiance or flux closely approximates the value of integrated irradi- 
ance over 0-90 degrees. The use of this angle then makes simple and 
inexpensive the computation of the net flux and consequently the atmospheric 
temperature change for a layer. The layer warming/cooling may be given by 

O (1440) , (10) 

where F is the net irradiance (cal cm'^min"!), p is atmospheric pressure 
(millibars), Cn is the specific heat at constant pnassure and 1440 is the 
number of minutes in a day. 

Water vapor inference . The following equations show the technique by 
which water vapor may be inferred from observed radiance. 

The modified equation (1) for downward radiance only is aiven by 

N+ = F / <1. (v) B(v,T(p)...) dpdv 

VI pr 

where u is the optical mass of water vapor. Here the downward radiance is 
measured, the filter function is known; a temperature profile above the 
reference is assumed by using the sounding of the nearest high altitude 
observing station. The altitude of the reference point (pressure) is 
known, and finally the frequency band of the radiometer is known, having 
been selected to be in the water vapor band. Consequently, the only unknown 
in Eq. (8) is the water vapor optical mass. 

Radiance may vary by changing the optical mass since 

T = x(w,k) (12) 

where k is the absorption coefficient for water. 
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An iterative method is then applied to equation (8) with a different 
value for u, until the difference between the observed radiance and the 
newly calculated radiance is minimized. This may be expressed as 

(N4-Q - < e (13) 

where e is the noise equivalent radiance of the radiometer. 

The optical mass varies with the reference level mixing ratio q^ in 
the mathematical approximation for u. This may be expressed as 

“ " 9 / ^ I %Pi *P- 

Pp r 1 

Here g is the acceleration of gravity, the subscript "r" refers to the 
reference level, and X is a power law exponent. By carefully choosing the 
value of X one may assume a constant mixing ratio or, as experimentally 
determined, a more realistic profile for the mixing ratio. 

P rofiling . The method by which the temperature or humidity profile 
from aircraft may be inferred utilizes the fact that the frequency of the 
filter or the angle of measurement (x) may be varied, thereby varying the 
level of the atmosphere from which the radiance measured originates. The 
upwelling radiance observed by the radiometer is a function of the atmos- 
pheric temperature and humidity profiles. 

The radiative transfer Eq. (1) may be rewritten as 



Here x refers to the frequency or angular characteristics of each 
observation^ (x,p^) is the atmospheric transmission between the aircraft 
pressure level (p^) and pressure p, and Pq is the surface pressure. B{x,p) 
is the Planck radiance and is given by 

B(x,p) = CivV[exp(C 2 v/T(p)-l] (16) 

where Ci and Cg are Planck constants. in (12) is called the 

weighting function. P 

Program RADIANCE uses a direct iterative method to obtain the inverse 
solution of (12) and (13), and the process is repeated until convergence 
between the observed and calculated radiance is achieved. This may be 
expressed as 

[!:(N(x,P^)caic - N(x,pp„|,5)2]^ < e (17) 

where e is the noise equivalent radiance of the radiometer. Rapid convergence 
is obtained by using a modified Newton-Raphson routine. 
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2.4 Infrared Radiation Observations and Calculations in the Atmosphere 
2.4.1 Background 

The SKYLAB field phase in June, August, and September Of 1973 in the 
Houston, White Sands, and Phoenix areas provided a unique opportunity to 
conduct infrared observations and subsequent calculations of the trans- 
mission and absorption properties of haze and middle altitude clouds. 

The observations used both an 8.0 to 13.0 ym downward looking chopper 
radiometer and an upward- downward looking radiometer (4.39 to 40 ym) 
attached to a standard radiosonde. 

One aspect of this research was to observe experimentally and analyze 
a volume absorption coefficient together with related optical properties of 
dust and middle tropospheric clouds to permit a simple solution for bulk 
radiation transfer phenomena. Since a detailed, moderately high resolution 
(1.0 cm“^ resolution) radiative transfer approximation program, RADIANCE, 
has been in use by the authors for several years, and since it being con- 
tinuously updated, it is only necessary to add these additional absorbers 
to the existing solution once their radiative character is determined. 

This in-situ method of determining the haze and cloud optical properties 
is different from other methods and therefore can provide interesting 
comparisons. 

For the radiometric transfer calculations, the simultaneous radi- 
ometersondes provided profiles of the free air temperature and humidity 
profiles. In each site the measurements were made over nearly homogeneous 
surfaces in order to provide a nearly uniform interface. 

The upwind and downwind locations chosen from some large source such 
as Phoenix enabled a determination of the IR transmission of both dust 
and altostratus clouds. In each case, the bases and tops were distinct. 


2.4.2 Some IR Radiative Properties of Haze and Clouds 

Platt (1973, 1974) as well as Kuhn et al. (1974) and Kuhn (1970) have 
investigated and reported their results concerning some of the optical 
properties of cirrus, contrails,, and middle altitude clouds. Preceding 
much of this work were the efforts of Hall (1968) in determining a physical 
model of radiative transfer in cirrus. The general agreement, at least 
insofar as the volume absorption coefficient is concerned, is very good. 

The investigation of middle altitude cloud radiative properties was prompted 
by a desire to compare these results with previously reported results of 
Platt (op. cit.) on middle latitude clouds. 

The monochromatic emissivity of a haze layer or a cloud layer on the 
basis of continuity is defined as 

° '• -"x - Px’ 
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where t is the transmission and p the IR reflectivity. From (1) we may 
write, ignoring p. 


Int^ = In(l-E^) = (-k^AZ), ( 19 ) 

where K (km“M is the volume absorption coefficient and AZ(km) is the cloud 
thickness or depth. We now define a bulk absorption coefficient, K^x> 
covering a small spatial interval, aa, from 9.5 to 11.5 ym hereafter symbol- 
ized as K(Kuhn et al . , 1975). 

To facilitate comparisons. Table 2 lists the haze and/or cloud physical 
properties with the symbols employed by the authors on the left and those 
of Platt (1974) on the right. 


Table 2. 'Physloal Pro-perties of Atmospherio Absorbers 


Th i s 

work 

Platt 


z 

h 

Optical Thickness : 

Au = / KAZ 

tA = / dh 


z=0 

h=0 

V 0 I. Absorp. Coef.: 

K(km"l) 

a(km"M 

Cloud Thickness : 

AZ(km) 

h(km) 

Thickness vs. e : 

1 

KAZ = In(l-e) 

= In(l-e) 


2.4.3 Methods of Analyses 

Consider a haze or cloud layer of thickness aZ, and optical thickness, 
KaZ, as in figure 1 . 


N2 



N22 
. — • — 22 


Figure 1. Symbolic atmospherio haze 
layers with arrows depicting up- 
welling IR radiance. Here t 
represents t., as above in (18). 


The radiance transmission (1-c) is given by 


^ N2+ -p(N?> -N^) 

N22'^ 

and the radiance emissivity is given by 

_ Njt -N2;i+ -p(N2+ -N^) 
^ ' Nt+ -N22+ 


( 20 ) 


( 21 ) 


where x, e, p, and N all represent parameters over the small wavelength 
interval, AX. Thus we treat a "bulk" x, p, e and K. The errors in the 
determination of the absorption coefficient are based on a radiometer noise 
equivalent radiance of 2.35xl0"®w cm~^ sr"^. This results in an average 
error of ± .05 in the transmission or emissivity, x and e, respectively. 

The error in the volume absorption coefficient from (2) is given by 


dK = 


1 ^ 

AZ X 


( 22 ) 


From Eq. (22) and from an average error of ±0.05 in transmission, the average 
error in the volume absorption coefficient, K, is ±.009 km“^. The flux or 
irradiance emissivity is required in place of the zenith emissivity for 
radiative flux transfer and IR cooling calculations. Elsasser (1960) shows 
that the flux emissivity, ef, is approximately equal to the zenith emissivity, 
e (our observations), of unidirectional radiance from a column of optical 
depth, KAZ, 1 .66 (KaZ). 

In these equations Nf is the upwelling radiance (w cm“^ sr"^), p the 
layer IR reflectivity, and the subscript "c" is the black body radiance at 
the mean temperature of the layer. Subscripts 1, 22, and 2 refer to layer 
levels (Fig. 1). "22" is the same level as "2" but is the upward radiance 

in a dust-free or cloudless atmosphere. For shallow layers Nif may be used 
in place of N 22 f in Eq. (20) with only a small error. However, in deeper 
layers {>1.5 km) N 22 t niust be employed. This necessitated a profile in a 
haze or dust-free area or outside of the cloud. 


2.4.4 Results 


Haze . Average values of the haze transmission calculated from observa- 
tions entered in Eq. (20), the volume absorption coefficient, K, and the 
optical depth, KaZ, are summarized in figure 2. These observations through 
Eq. (19) provide the curve of K. The optical thickness, KaZ, is calculated 
and appears at integer values of the geometric depth of the dust. The mean 
value of K is 0.042 ± 0.009 km~^. However, the rms deviation is only 
0.0035 km-i. 

The smooth nature of the transmission and more important of the ab- 
sorption coefficient suggests a uniform distribution of the atmospheric 
dust load. Error limits of observations and calculations are given’by the 
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Haze Depth (h m) 


F'igio'e 2. Haze transm-issiorii volume absorption aoeffiaient relative to haze 
depth and optical depth. Data used ove from White Sands ^ Weu) Mexico y 11 and 
12 August 1973. 


the vertical arrows. The narrow range of variations in these optical prop- 
erties is in contrast to the much wider range of variation in clouds. Haze 
IR properties are an order of magnitude less, numerically, than those of 
middle clouds. The assumption is further made that the dust particles are 
uniformly less than a micron in diameter. 

Middle Clouds . Table 3 lists the observed and derived optical prop- 
erties of tropical middle altitude clouds and may be compared with the prior 
research of Platt (1974). In fact the cloud observations and calculations 
were made as a result of Platt's excellent work with hiS lidar. 

In the last two sets of As values, the large values of K, the high 
emission, and the strong absorption indicate that the clouds were composed 
of large numbers of supercooled water droplets. It is also evident that a 
mixture of ice crystals and water as well as various particle densities 
causes the range in the values of K. However, there is a tendency to relate 
moderately deep, ice crystal clouds with K values <1.0 km"^ and shallow 
dense altostratus, water droplet clouds with K values >7. 0 km-^. They are 
apparently composed mainly of ice particles with some supercooled water drop- 
lets evidenced by their lower values of K(K<T. 00 km) as suggested by Platt 
(op. cit. ) . 
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Table 3. Observed and Derived Cloud Optieal Properties 



Depth 

(km) 

Mean 

Mean 




Cloud Type 

Height 

' Temp 
(°C); 

T* 


K 

As-Ac 

1.5 

6.4 

-n.8 

0.42 

0.58 


As 

0.6 

6.25 

-11.0 

0.63 

0.37 

0.77 

As 

2.0 

6.4 

-11.8 

0.25 

0.75 

0.69 

As-Ac 

2.0 

6.2 

-10.0 


0.77 

0.73 

As 

0.4 

5.7 

- 7.0 

0.04 

0.96 

8.06 

As 

0.3 

5.2 

- 4.5 

0.07 

0.93 

8.86 

*Cal culated 

from Eq. (3) 

as an average 

of ten 

observations 

of N2f, 

N22f and 


N +. 
c 


2.4.5 Calculated Cooling Rates 

Atmospheric Haze . The NOAA RADIANCE program discussed in Section 2.3 
has provisions for computations involving various aerosol layers such as clmH 
or haze, providing the volume absorption coefficient is known. 

To calculate cooling rates through a haze layer, a temperature/pressure/ 
humidity profile over White Sands, New Mexico, for 11 August 1973 was input 
to the RADIANCE solution. Figure 3 is a plot of the results, which provide 
calculated IR cooling rates at the midpoints of each layer in the haze and 
a corresponding haze-free and cloudless cooling rate profile. 

The computations of infrared cooling show the cooling to average only 
0.009C hr"i throughout the haze layer. An absorption coefficient of 
0.042 km-^, as determined from observations, was used (Fig. 2). This cooling 
may be easily balanced by solar warming in the layer. This cooling rate 
would hardly alter the dynamics of the haze layer. 

Middle Clouds . Radiometersonde profiles through middle altitude 
altostratus and altocumulus clouds furnished the RADIANCE calculation input 
to determine cooling rates through these clouds. Absorption coefficients 
from Table 2 were employed in two examples with calculated cooling rate curves 
appearing in figure 4. The curves are labeled, cloudless,; K=0.77 and K=0.8. 
Cloud bases and tops are indicated for the two calculations. 

Moderately deep clouds with K<1.00 exhibit calculated cool ing rates of 
up to 0.23C hr"^, not enough to alter the cloud stability. In contrast, 
shallow dense clouds of water droplets with K = 8.0 display calculated 
cooling rates of approximately 0.15C hr“i. Such a cooling rate could have 
an effect on cloud stability. 
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-O.'l -0.2 


'C/hr 

0 


0,2 0,4 



F'lgupe 3. Cdtoulated IR aoot-ing in Figure 4. Calculated IR cooling 

haze^ 21 August 1973, 2130 UT, through thin cdnd deep middle clouds, 

over sand dunes, Iflhite Sands, 8 August 1973, 1937 VT over 

New Mesd-co, Rosenburg, Texas, oil fields. 


2.4.6 Conclusion 

It would appear important now to combine IR cooling rates for atmospheric 
absorbers with solar heating to complete the radiation budget analyses for 
both cloudless and cloudy and hazy atmospheres. Also the bulk radiative 
transfer approach employing a volume absorption coefficient determined in 
the atmosphere for non-gaseous absorbers and combined with a mass absorption 
coefficient (cm^g'i) for the gaseous absorbers provides a simple and rapid 
computer solution. In cloud or in haze, scatter of radiant emission was 
assumed negligible but the nature of the observations would appear to validate 
the use of the experimentally determined absorption coefficient in the 
presence of haze or cloud. A classification of absorption coefficient in 
relation to various cloud types and depths as well as to various haze or 
dust layers could now be attempted from the large amount of data available. 
Further, sufficient simultaneous haze and cloud composition and IR observa- 
tions exist to attempt this classification as input to radiative transfer 
calculations for future improved input to numerical weather models. 
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0 




Figure 5. Profiles of surfaoe temper- Figure 6. Caloulated brightness 
ature from helicopter airborne. terrperatures . 

observations 3 and as calculated with 
a transfer model , using Flsasser 's 
table (K^)j Smith's table (^^(s))^ 
and this study 's Phoenix bulk co- 
efficient (Kpyx)> 


The special effects of a heavily dust laden atmosphere on the magnitude 
of the bulk total extinction coefficient for the 11.0 - 11.1 ym spectral 
interval are evident after a review of figure 5. This figure illustrates 
calculated and observed brightness (apparent) surface temperatures as a 
? function of height over Rainbow Valley (Phoenix, Arizona). The differences 
f among the calculated and observed profiles clearly demonstrate the necessity 
- of obtaining "K" for a unique area. 

Figure 6 illustrates the calculated brightness temperatures at 10.0 and 
12.5 ym and at 8 to 14 ym. A computation using 11.0 to 11.1 ym gives a 
curve close to the satellite curve. It also agrees with calculations made 
with the Phoenix area determined total extinction coefficient. 

Figures 7 and 8 show the sites of the measurements used in figures 5 
and 6. Figure 9 summarizes the various bulk total extinction coefficients 
employed in calculations. Kjg-j- for this paper provides the best results when 
compared with S-191, 11.0 - iT.l ym channel. Note that v (cm* ^) is the 
inverse of the wavelength in microns. 
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SKYUI MAXIMUM ZOOM ^HOTOMA^ 
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Figure 7. SKYLAB photomap of Phoenix ^ Arizona ^ area 
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Figure 9. Total extinction coefficients (cm^g~^) according to Bignell (K^)j 
Kuhn utilizing Smith CO2 tables (K(S)), Kuhn utilizing Elsasser CO 2 tables (E)j 
and this study (PHX). 


2.5 Absorption of Tropospheric Aerosols: Urban and Rural Aerosols 

of Phoenix, Arizona 

The infrared optical properties of aerosols around Phoenix, Arizona, 
were studied to evaluate the effects of the interaction of atmospheric 
particulate matter with long wave radiation. This was accomplished by 
collecting aerosols on membrane filters at 500 m msl that had air asperated 
through them at a controlled rate. A Gardener fine particle counter was 
run concurrently. Filter samples were then analyzed for particle sizes by 
a scanning electronic microscope. The aerosoT absorption coefficient was' 
calculated according to Mie theory. By measuring the aerosol layer trans- 
mission, the bulk total extinction coefficient was determined. Further 
calculations show that at 10 ym wavelength, more than 95 percent of the 
infrared extinction is caused by absorption. The imaginary part of the 
infrared refractive index of 0.47 for urban and 0.19 for rural aerosols was 
calculated by Pueschel and Kuhn. (1975). V 
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2.6 Discussion and Conclusions of the Surface and Boundary Layer Experiment 

This field program, under EREP SKYLAB overpasses of selected sites, was 
designed to achieve three primary objectives: 

1. Provide in-situ balloon-borne radiometersonde observations of . 
total (infrared plus solar) upward and downward irrariiance (flux) 
from the surface to 35 km for comparison with and verifications of 
selected atmospheric radiative transfer models. These results 
were correlated with the S-191 broad-band data. 

2. Determine the mass absorption coefficient (cm^q”^) in the atmos- 
pheric window (8.0 to 14.0 ym, 714 to 1250 cm"^) by direct 
relating of S-191 11.0 to 11,1 ym (900 to 910 cm“^) brightness 
temperature observations to surface and helicopter radiometric 
observations in the same spectral band. 

3. Provide interface (air-surface) observations of net total radiation, 
net total insolation reflected insolation, downwar<i insolation 

and albedo for calibration, ground baseline, and comnarison with 
S-191, S-192, airborne, and balloon-borne radiometric observations. 

Conclusions relative to the first objective were that the balloon 
radiometersonde, a device flown successfully since 1957, does provide a 
comparison for radiative transfer models that is at least as accurate as 
the models. Deviations of radiative flux and intensity calculations are 
due to inadequate knowledge of real atmosphere transmission. The SKYLAB 
S-191 radiometer was vital in providing correlative atmospheric long wave 
intensities (radiance) in the center of the atmospheric window region. 

Relative to, the second objective, it was possible to determine the 
mass absorption coefficient of the atmosphere over at least the Phoenix site 
in the 8 to 14.0 ym (714 to 1250 cm- ^) spectral band and in the 11.0 to 
11.1 ym (900 to 910 cm"^) spectral band by combining observed S-191 bright- 
ness (radiance) temperatures with balloon radiometric and helicopter 
radiometer data. This results from in-situ field observations requiring 
the S-191 observations is important in that it provides parameters required 
for proper atmospheric radiative transfer under a typical desert southwest 
atmosphere under summer conditions. 

The third objective, that of providing ground baseline observations 
for S-191 and S-192, basically a data-base effort, was also achieved. 


3. AIRCRAFT EXPERIMENT 

3.1 Comparison of Models For Computing Atmospheric Infrared Transmission 

Various models were used in describing the primary differences in 
transmissiveness of gaseous constituents of the atmosphere. Among them 
were RADIANV of Colorado State University and RADIANCE of the NOAA, 
Environmental Research Laboratories. These include options for computing 
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radaiation transmission including ozone, aerosol, and continuum as well 
as water vapor and carbon dioxide transmissions. 

These models were then compared with SKYLAB infrared spectrometer 
data. It was shown that all models are capable of computing atmosnheric 
absorption and emission from the earth's surface to within 5 percent 
accuracy. The models are described in detail in Renne and Marlatt, 1975. 

3.2 Spectral Estimates of Albedo and Comparison to SKYLAB Observations 

An application of the "adding" or "doubling" method has been made for 
homogeneous atmospheric layers composing an inhomogeneous atmosphere 
(Rainey and Marlatt, 1975). The lower layer was assumed to be composed 
of silicate particles of complex refraction index 1.55 - Oi corresponding 
to a wave length of 0.5 ym. The particle size distribution was measured 
on site. The upper layer was assumed to be a gaseous layer scattering 
light in accordance with Rayleigh's phase function. Radiance Values were 
computed for a planetary system composed of the described atmospheric 
layer overlying a Lambert reflecting ground. These values show a fairly 
good agreement with the SKYLAB observations given a reasonable knowledge 
of the scattering (or reflecting) properties of the atmosphere and ground. 

3.3 Conclusions 

This field program at aircraft level was designed to determine the 
accuracy and applicability of a number of radiative transfer models for the 
atmosphere by comparing. values from the models in conjunction with observa- 
tions from the SKYLAB spacecraft, from balloon radiometersondes, and from 
surface radiometers. 

Note that the transfer models studied covered, basically, the same 
frequencies of the spectrum as were employed in the S-191 and S-192 SKYLAB 
radiometers. In fact the 11.0 to 11.1 ym (900 - 910 cm"^) band of S-191 
was basic to transfer calculations in this atmospheric window region, 
providing interface brightness (apparent) and physical temperatures. The 
S-192 data provided the areal background imagery for verification. 

The results of this study were: (1) Each of the models tested was 

able to compute the long wave radiation emitted from the earth's surface 
to an accuracy of 5 percent of observed radiance. (2) The models were 
all most sensitive to input parameters, especially surface temperature, 
water vapor and ozone. (3) Particularly in dry atmospheres, all models 
tended to overestimate the opacity of the atmosphere in the 8-14 ym band- 
pass. 

4. DISCUSSION AND CONCLUSIONS DERIVED FROM SCARP 

The SKYLAB Concentrated Atmospheric Radiation Project succeeded in 
utilizing radiometric observations from the SKYLAB Earth Resources Experi- 
ment Package. Its primary goal was to arrive at a more complete description 
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and model of radiation transfer within the atmosphere by employing in-situ 
satellite, aircraft, ground, and balloon-borne radiometric observations. 
Specifically an attempt was made to determine the applicability of a variety 
of radiative transfer calculation models computing radiative flux in air 
masses of different character, i.e., varying amounts of wet and dry 
aerosol burden, varying interface types, and varying temperature and 
humidity profiles. 

The method of comparison was one of obtaining in-situ SKYLAB S-191 and 
S-192 observations, aircraft and balloon-borne profiles of upwelling and 
downwelling radiance and irradiance, and ground-base observations of net and 
total infrared and solar radiation. Having temperature and humidity 
sounding profiles under the EREP overpasses, we ran the transfer calcula- 
tions in both the solar and infrared spectra and compared them with the 
in-situ observations. The field data for comparisons were good and these 
data sets are available for furthei^ research, either as computer plotted 
graphs in this report or as card ard/or tape data sets. Printer output is 
also readily available. From this basic effort one may conclude that, 
direct observations are superior to present transfer calcinations, and 
this is because of a lack of knowledge of complex atmospheric transmission. 

This comparison of several existing infrared radiative transfer models 
under somewhat controlled conditions and with the above atmosphere ob- 
servations of SKYLAB's S-191 and S-192 radiometers illustrated that the 
models tend to over-compute atmospheric attentuation in the window region 
of the atmospheric infrared spectra. Add to this the variability of radiance 
calculations from the different models, all with identical input data, and 
one is led to the conclusion that much further research on in-situ atmos- 
pheric tv’ansmission observations and calculations is a necessity. This is 
further evidenced by inaccuracies in satellite temperature and humidity 
profile specifications from an inverse solution to the transfer equation. 

In this connection atmospheric transmission functions for at least the 
Phoenix summer atmosphere appear to have been determined. 

A basic recommendation for future studies would be to employ the SCARP 
accumulated satellite, aircraft, and balloon platform acquired data in an 
attempt to determine atmospheric, broad-band transmission coefficients for 
the Phoenix and Houston areas. This could be accomplished from the data 
set on hand and could serve the radiative transfer modeling community well. 
Such an undertaking could be undertaken for little more than computer costs 
and some salary budgeting. 
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Plots of Observed Radiation Flux Profiles 
During SKYLAB 
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APPENDIX B 


Data Listings of Observed Radiation Profiles 
During SKYLAB 
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U. 3i; 47 

-8.0 

A 


0. 38.1 

0 . 0657 

0. 3144 

-8.3 

A 


0.3958 

0. .j6 31 

0. 32 67 

-5.2 

A 


a.39u5 

C. 06 71 

0. 3294 

-1.1 

A 

*■ 

0. 38 76 

C. 06C9 

0. 3267 

2.3 

A 


a .-3 7 68 

0.0557 

0.3212 

3.0 

A 


0. 37 25 

0. 0567 

0. 3158 

-0.8 

A 


0. 3805 

0.0624 

C. 3180 

-9.6 

A 


a. 39 51 

0.0637 

0.3314 

-16.8 

A 


0. 4012 

0 . 0633 

0. 3378 

-10 .9 

A 


0.3996 

0.0623 

0. 5373 

-5.3 

A 


0. 40 30 

0.0628 

Q.34G2 

-5.3 

A 


0.41U2 

0.06^3 

0. 3459 

-7.4 

A 



,5 

0. 41 63 

0.0636 

0. 3521 

-9 

>5 

0. 42 3-. 

0 . 0 b 64 

Q. 3570 

-9 

,0 

G. 4345 

0. 0715 

0.36 30 

-7 

,4 

0. 4399 

0 . 0750 

C • 3d49 

-7 

.5 

0.44C5 

C. 0761 

0.3645 

-8 

,4 

U. 4427 

C, 0685 

0.3741 

-7 

,7 

0. 44 35 

C. 0654 

0.3781 

-4 

.0 

0. 44 86 

0. 0709 

0.3777 

-0 

.9 

0.4529 

G. J786 

0. 3744 

-1- 


0.-5 74 

0.0751 

C.37B3 

-5 

.3 

0. t-5 91 

0., 07 67 

0.5823 

-6 

C 

0.-649 

C • Ofl G? 

0. 3843 



























3^2.2 
3 30 .3 
318.1 
305. 5 
2 92.3 


280 .3 
269.7 
2 58.8 
2U7.3 
’235.8 


2 15 . 3 
2 06.9 


189.7, 

130.3 

172.2 


164.3 

156.8 

149.5 


0 . f 1 46 
’ .0.7131 


0. 70 83 
-24.7 Q.7Q89 

- 2‘6 .7^~T,7l^5 
-28.9 _J.71C3 
- 31.2 0, 70 67 

-33,4 0.70 38 


0. 5797 
C.5791 


0. 5738 
JO. 5647 
o'. 55 59 
0.5400 
0. 5297 
0 .53 05 


C. 13 ‘.a 
0.1540 


0.1350 
C. 1h42_ 
d'7l5 46 
0. 1703 
0. 1770*“ 
0. 17 34 


2.0 99 
1.841 


1.556 

1.259 


88 
0 . 861 


0.6 
J.5 32 
0.445*' 
0.373 
*0.277 
0.236 




0. 7067 

0. 52CC 

0.18 67 

-2.4 

0.199 

30 

0, 7074 

Q.51ii 

G, 1963 

-0.5 

0.159 

31 

0. 7065 

0.5177 

0.1838 

2.5 

0.127 

30 

0. 70 50 

0. 5275 

0. 1775 

3.6 

0.102 

30 

0.7083 

0. 5359 

0.1724 

2.1 

0.0 85 

31 

0.7131 

C . 5435 

0.1696 

-1.9 

0.Q65 

32 


163 

0 . 29 

0.17 35 

-7.7 

A 

199 

0.52 70 

0. 19 30 

-10.3 

A 

2 29 

0.50 17 

0. 2212 

-6.2 

A 

2‘+0 

0.5027 

0.2213 

-4.3 

A 

2 23 

0.5024 

0. 2199 

-3.4 

A 

2 32 

U « ( G 

0.22bl 

-5.1 

A 


0 

-bC ,8 

0 

-62.3 

0 

- 6i* . 4 
















IFILTERED 

kAljlOMcl 

RUN FOP 

RN9 

PRE SS 

TT''*5 

t-A IK 

r-UP 

F-OM 

(HP,) 

(MIN> 

(IFG GMLY/HIN) (LY/HINJ 

lOlt.7 

• 0.0 

33 .1 

0. .7173 

0. 7011 

9-^9. 3 

1.0 

25.0 

0. 80 31 

b. 7214 

963.1 

2.C 

21 .0 

0.7816 

C.7179 

937.3 

3.0 

19.4 

0.7674 

0.7178 

910 .8 

4.0 

17 .7 

0.7675 

C i 71 63 

8 66.6 

^ ■ G 

lb. 2 

0.7698 

u . 7 0 0 5 

8 63.3 

6.0 

15 .2 

0. 7785 

C . 67 96 

6 40 . 3 

7.0 

14.2 

0. 7840 

0. 6742 

817.3 

^^KKI 

12.9 

0.73 33 

0.6771 

7 94.5 

9.0 

11 .6 

0.7749 

0 . 6738 

771 .3 

10.0 

10.1 

0.76 33 

0i6679 


11.0 

6.6 

0.75 95 

0.6501 

7 26.1 

12. u 

7.2 


0.6432 

704.7 

13.0 

5 .8 

0.7495 

0 ■ 64 Z 3 

6 84.3 

l4.C 

4.7 

0.7414 

C . 63 36 

664.8 

■BSQ 

3.3 

0. 73 64 

0.6311 

646.3 

16 . 0 

2.0 

0.7279 

0.6290 

6 26.5 

17.0 

. C .7 

0.7227 

0.6185 

605.3 

18.0 

-0 .4 

0.7176 

0.60 96 

583.3 

■jnsHM 

-1.9 

u. 71 67 

0.6053 

5 64.2 

20.0 

-3.3 

0. 70 92 

0.6020 

548 .2 

21.0 

-4.7 

0.7056 

0. 5949 

5 32 .3 

22.0 

-5.6 

0.7006 

0.5855 

515.5 

23.0 

-6.4 

0. 6973 

0.5760 

497.6 

24.0 

-7.6 

0 . 69 26 

0.5711 

4 807? 

25.0 

-9 .3 

0.6883 

0.5669 

463.7 

26.0 

-10 .9 

C . 63 64 

0.5606 

4 47.0 

27.0 

-12.6 

C. 6 846 

C . 5549 

429.7 


-14.5 

0.6836 

0 . 5522 

415.1 

29.0 

-16 .1 

0.6811 

0 . 54 71 

403.2 


-17.7 

0. 67.91 

0.5363 

392.4 

31.0 

-19.2 

0.6774 

U.5261 

378.2 

32.0 

-20.6 

0. 6748 

0.5219 

361 . G 

■■ 33.0 

-22.3 

0.6689 

'0. 5190" 

343 .6 

34.0 

-24.5 

0.6681 

0. 5138 

330 .3 

35.0" 

4 26.4 

0 . 66 74 

C.5Q88 

319.2 

36.0 

-28 .3 

0 • 66 46 

0 . 49 84 

307.1 


-30 .2 

0.6662 

0. 4370 

294.1 

^khII 

-32.5 

0.6725 

0.4878 

233 .5 

39. C 

-34 .8 

0.6795 

0.4940 

274.7 

40.0 

-36.3 

0.68C5 

0.4902 

2 64 . j 

41.0 

-39.3 

0. 6313 

0. 4859 

252.3 

42.0 

-41.7 

0.6904 

C.4899 

HKCIflHSi 

43.0 

-45 ,9 

3. 70-35 

0 . 49d4 

229.7 

-4.0 

-45.7 

U. 72 73 

0 . 49 69 

2 20 . 1 

45.0 

--7.5 

0. 7383 

0. 4900 

211.3 

,46.0 

-49 .6 

Q. 7391 

G • 48 3Z 

202.6 

47.0 

-51.7 

0. 75 40 

0.4773 

193.1 

48.0 

-53.5 

0. 7661 

0. 4774 

183.1 

. 49.0 

-54.9 

0.7357 

0.4779 

174.0 

50. Q 

-56.8 

0. 30 86 

0. 47 84 

166.2 

51.0 

-5 8.6 

"Q.'ass?' 

0.4814 

158.7 

52. 0 

-5P.7 

0. 85 77 

C. 4810 

150 .9 

53.0 ' 


0. 8813 

0.4777 

143.4 

54.0 

-62.7 

0. 90 69 

0. -779 

136 .6 

'55 .0 

“ 04 . 3 

0.9398 

C. 48 24 

130.2 

56.!: 

- 64 . 6 

0. 9653 

0. 4363 

124.3 

5 7.0 

-63 .9 

Q-, 9?22 

0 . 48 eo 

118.5 

5 8.0 

- 62 .5 

0.97=0 

C. -867 

112.3 

59.0 

-61.4 

0*9791 

0. 4856 


60.0 

-.60 .8 

0.98 53 .■ 

0.48C6 
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3 AUG 

1973 

OPUS 11 

F-NCI 

COOL 

Q-MIX 

R-HU'I 

(LY/MIN) 

<0£G/nA) 

(G'!/KG) 

(PC) 

0.Ulr2 

0 .0 

23.590 

83 

"O.'OHi? 

-I5i0 

14.524 

76 

0.06 37 

1.7 

14.0 O', 

86 

0. 0 4 9b" 

0.9 

13.305' 

86 

0.0512 

-2.1 

10 .7 12 

77 

0.0 6 93* 

-4‘. 1 

9.704 

75 

0.0983 

-3.3 

9.737 

76 


-1.8 

8 . 0 31 

67 

0.1062 



70 

""o.ioii 


7.643 

70 


-0,3 

7.195 

71 

0.1094 

-0.4 

"6.719 

71 

0,1068 

-0.0 

6.5 36 

74 

0 . 1 G 69 


6. 0 21 

73 

0.1078 

0.6 

5.772 

74 

OTlO 53 

F. 4 

5.F4Q 

79 

0.0989 

.0.0 

5.3 69 

80 

o'. id-£ 


5'.426 

82 

0.1080 

-0.7 


86 

0.1114 

-0.4 


76 

0.1072 

-0.4 

3. 845 

72 

C. 1107 

-f.3 


72 

0*1151 

-1.4 

3.260 

68 

0 .1218 


2.8 76 

61 

0.1216 

-0.7 

2.5 53 

58 

0.121- 

-0.7 

2.143 

54 

0.1257 


1.768 

49 

0.12'98 

-1.1 

1.T37 

53 

0.1313 

-1.5 

1.6 30 

'57 

o'.Tsvb 

- 2 .“4 

1.413 

54 

0.1429 

-2.9 

0.987 

41 

0.15 13 

-1.9 

0. 983 

46 

0.15 29 

-0,8 

0.778 

39 

0.1499' 

“ ■-0 .6 

0V623 

35 

0.1543 

-1.4 

0.523 

34 

'■“0,15 85 “ 

-3.2 

0.442' 

53 

0,1663 


0 .318 

26 

8*1792 

-3.6 


26 

. 18 46 

-2.8 


26 

0.18 56 

-2.F 


- 24" 

0.1903 

-2.4 


25 

0.1959 

-3 . 4 

0.134 

26 

0. 2CC5 

-4.9 

0.100 

27 

5. 2122 

-7.1 


25 

0.2304 

-8.5 

0..0 67 

24 

0. 2487 

-9". 3 

0.0 55 

24 

0. 25 58 

-9.1 

0.0 55 

26 

~*0'. 27 27" 

-9.6 


A 

0. 2887 

-11.5 


A 

0.3078 



A 

__0. 3 3 02 

-15.4 


A 

0,35-3 

■"-17.7 


a' 

0. 37 67 

-18.8 


A 

■"0."4b 36 

-20.5 


A 

0 ; 42 90 

-21.1 


A 

0.4573 

-16. 8 


A 

0. 4790 

-14.1 


A 

0 . 43 55 

-7.9 


A 

0.4683 • 

-5.8 


A 

0* 49 35”“ 

-4.7 


A 


0. 5047 



















0. 9816 

0. h799 

0.5013 


0.1,356 

G . -.9 61 

■3 . 97 35 0 . £. 8 95 

0,^395 

C,.%9,9 

. C.* o » ' 

Q. Vff53 


S .3 ■ d. «550 

0. ^7:94 

l).'^7 56 

6.7 0.939*+ 

0.47 61 

0. 46 33 

6.5 0.9'2 8S' 

0. U773 

0. U6U7 

6.7 Q.9262 

C .43*,Q 

0.44 22 


a. 9911 

0. 48 96 

G.'i+GiS 

0.8863 

C. 4928 

0.39 35 

0. 8828 

■Q. 49 56 

0. 3971 

a. 96 76 

0.4944 

0. 37 32 

0. 8514 

Q. 4924 

0.3589 

0.8402 

G . 49 7C 

0 • ^ o2 


0. U3 64 


0.79 13 

0.5 a 30 

0.2883 

~u, 7914 

c. sorg" 

■"0.~28'35“ 

0^7 9 j) 2 _ 

0.5093 

0.2809 

0 . 78 56 

0 • 50 68 

0.27 87 

0.7785 

0.5024 

0. 2761 


i 











IFIL-HKEO 

4ADI0METERS0NDE 

RUN FOR 

RN3 ■ 

3 AUG 

1973 

OPUS 

12 

PRESS 

TIPE 

T-AIR 

F-UP 

F-ON 

F-NET 

COOL 

O-HIX 

R- 

HUH 


( MItl) 

(DEG C) (LY/M.IN) (LY/riK) (LY/HIN) 

(OEG/OA) 

(GM/KG) 

(PC) 

10 13.1 

0 . 3 

31.5 

0.7024 

0. 6904 

0.0120 

0.0 

22.361 


69 

99C.6 

1.0 

27 • 6 

G. 8470 

0.6489 

0.1961 

-48.6 

17.194 


72 

96*r.9 

2.0 

25.6 

0. 6343 

0. 6371 

0.1972 

-0.9 

15.755 


73 

937.2 

3.0 

23.1 

0. 6167 

0. 6275 

0.1912 

1.2 

14.364 


76 

904.5 

4.0 

20 .6 

0. 80 45 

0.624J 

0.1805 

1.0 

15.056 


87 

375.1 

5.0 

‘ 19 .3 

0 . 7946 

0.6160 

0.1766 

0.4 

14.712 


92 

6-»9.5 

6.0 

17.6 

0. 7874 

0.6094 

0.1780 

-C .3 

11.695 


76 

622.9 

7.0 

15,6 

0.7803 

0 . 5964 

0.1836 

-0.8 

10.976 


81 

795.3 

3.C 

13.7 

0. 7811 

0 . 5968 

0. 18t3 

-0.7 

11.501 


94 

767.7 

a. 3 

12.4 

0.7324 

0.5915 

0.1969 

-G.l 

11.059 


92 

. 740.0 

10.0 

11.1 

a. 7778 

0. 5367 

0.1911 

0.8 

9.282 


83 

713.3 

11.0 

9.2 

0. 7656 

0. 5838 

0.1818 

1.8 

8.658 


64 

669,0 

12.0 

7.2 

0 . 7 6 35 

0. 5938 

0. 1696 

1.8 

7.754 


83 

667.0 

13.0 

5.4 

0.76«2 

0. 6056 

0.16 34 

0.1 

6.909 


83 

645.6 

14.0 

4.1 

0.7814 

D. 6170 

0.1645 

-3.7 

6.783 


84 

625.1 

15.0 

2.9 

0. e023 

0.6205 

0. 1821 

-8.1 

7. 006 


92 

605.3 

16.0 

1 .4 

0.8373 

0. 6108 

0.22 62 

-10.6 

6.190 


89 

535.5 

17.0 

C .2 

0. H678 

0.5959 

0. 2720 

-10 .5 

4.745 


71 

5 65 .3 

16.0 

-1.2 

0. 8844 

C. 5641 

0.3002 

-8.0 

4.510 


72 

545.2 

19.0 

-2.6 

0. 8958 

0.5727 

0. 32 32 

-1,5 

4.297 


75 

525.7 

2 u . 0 

-3.9 

0.90 95 

0.5743 

0.3351 

6.4 

4.050 


72 

5 06.1 

21.0 

-5 .4 

0.90 84 

0. 6279 

C, 2805 

18.2 

3.606 


72 

437 .3 

22.0 

-7 ,4 

0.8945 

0.7123 

0.18 22 

20.2 

3.006 


65 

469.7 

23.0 

-9.Q 

0. 8626 

0. 7757 

0 . 1 0 69 

12.5 

2 . 4 76 


61 

453.0 

24. C 

-1C .7 

0. 8826 

0. 7740 

Cl . 10 88 

1.1 

1.994 


52 

437.3 

25 .C 

-12.5 

0. 8933 

0. 75 93 

0.1350 

-5.7 

1.789 


54 

422.7 

26.0 

-14.5 

0. 90 81 

Q . 7544 

0.1536 

-4.0 

1.595 


54 

408 .2 

27.0 

-16 .5 

C. 9253 

0 • 7664 

0. 1586 

0.4 

1. 412 


53 

3 92.9 

28.0 

-16 .3 

0.9341 

0.7864 

0.1477 

Z 90 

1.122 


50 

376.9 

29.0 

-20 .3 

0.93c5 

C. 60 32 

0.1333 

-0.8 

0.937 


45 

360.1 

80 . 0 

-22.2 

0.93 28 

G . 8 0 j 7 

0. 1322 

-11. 1 

0 . 656 


37 

342.7 

31.0 

-2v.8 

0. 9130 

0. 7‘.13 

0.1767 

-25.6 

0.423 


28 

3 26.3 

52.0 

-27.4 

0.9139 

0. 6293 

0. 28 40 

-37.1 

0.422 


34 

312.1 

33.0 

-29.7 

0.9172 

U.5031 

0. 4140 

-37.7 

0.313 


30 

299.8 

34.0 

-32.1 

0.9247 

C, 4342 

0.4905 

-26.8 

0. 214 


24 

2 66.0 

35.0 

-34.7 

0.91c3 

0.4064 

0.5099 

-12.6 

0 .228 


32 

2 75,7 

36.0 

-37 .2 

0. 91u2 

0.3910 

0. 5192 

"6 . 1 

0.231 


40 

2 64.7 

37.0 

-39.2 

C.9046 

0. 3762 

0.52 84 

-8.1 

0.213 


44 

255 .7 

3 8.0 

-•41.0 

0.90-.9 

0,3652 

0 . 5 3 97 

-14.0 

0.199 


47 

2 47.9 

39.0 

-42.7 

0.92 54 

Q » 3 ^ 8 d 

0 . 5 6 86 

-23.8 

0.183 


50 

239.2 

40.0 

-44.9 

0,9664 

0 . 3598 

0,6065 

-33.6 

0.153 


51 

228.4 

41.0 

-47 .6 

1.04 47 

0. 3733 

0.6714 

-46.0 

0.117 


50 

216,6 

42.0 

-50.9 

1.1656 

0.3946 

0.77 10 

-57.4 


A 


2 04.7 

43.0 

- 5 3 . 5 

1. 32 82 

0. 4193 

0. 9129 

-59.5 


A 


194.3 

4^.0 

-56.0 

1. 4754 

0. 4330 

1. 04 25 

-41.5 


A 


164 .7 

45.0 

-58.6 

1. 5503 

0.4584 

1. 09 19 

-5.2 


A 


176.0 

HO 

- 6 C . 6 

1. 5261 

a . 49 1 3 

1 . u 0 43 

25.3 


A 


166.3 

47.0 

-61.8 

1. t795 

0.5220 

0,95 76 

31.0 


A 


161.7 

46.0 

-62.8 

1. 4769 

0.5417 

0.9352 

15.9 


A 


155.1 

4 9 . C 

- 64 . 9 

1. 50C1 

C. 55 36 

0.9465 

4.8 


A 


1 47 . 9 

50.0 

-66.5 

1. 50 30 

0.5611 

0 , 94 o3 

5.0 

* 

A 


141.3 

51.0 

-67.3 

1. 4859 

0.5619 

0. 9240 

8.1 


A 


154,5 

52.0 

-67.7 

1. 4780 

0 . 5606 

0.9175 

5.4 


A 

* 

127.2 

5 3.0 

-67.7 

1. 47 39 

D.56Cil 

0.9133 

1.0 


A 


120.3 

5 u . Q 

-67.6 

1. 47 8G 

0.5584 

G.9196 

2.3 


A 


114,9 

6 5.0 

“67,6 

1, 47 67 

0 . 5 5 9 3 

0.9174 

9.7 


A 


109.9 

56.0 

-67.9 

1. 4615 

0. 5597 

0.9016 

16.7 


A 


10 4.0 

5 7.0 

- 6 3 .0 

1.4399 

0.5632 

0. 87o7 

13.9 


A 


98.2 

56.0 

- 67 . 6 

1.4251 

0.5632 

G. 8619 

8.3 


A 


93.2 

5 9.0“ 

-67,4 

1. 4322 

0.5603 

T78719 

4.3 


A 


99.0 

60.0 

-66.0 

1. 4228 

0. 5545 

0. 66 82 

5.2 


A 



51 


I . > 


■ 84.9 

61.0 

••63 .6 

'1.4 063 

O'. 549d 

0 . 8563- 

5.8 

_ _ ^ 

80 .7 

62.0 

“61.3 

1. 3972 

C . 5464 

0.6508 

0.3 

A 

■” 76.3 - 

63 . j 

” 0 0 .8 

'1. 4Q 23 

0.-j4 2d 

0 . 8597 

2.2 

A 

71.8 

61.. 0 

“01.3 

.1. 3998 

0. ?360 
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-.3 

16.6 
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56 
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-2.6 
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bb 

74 
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“42~ 
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-1.1 

10.417 

87 

“ rj.“5“~0.o478 rri9iy~o.'o54i. 
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11.8 u.5463 0. 
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5497 0.C440 
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-214 0.1354; 

4066 0.14i4 

3928 071595 

3842 0.1657 


• 0 . 1 _ 

“5.0 7.538' 

0.0 6.273 

0.9 5.676 

2.9 5.253 




4.919 

4.328“ 

3.880 

“3.-55 
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•'19 .“3 07539J 0.3629 0.l7t4 
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Subroutine TABPRINT 

Entry Points: VLNPRNT, TWATER, & TCARBON 
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Subroutine PTPLOT 

Entry Points: UNIT & TPLOT 


Subroutine SORT 
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Subroutine HEADWAT 
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Subroutine CALWAT 
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o 0]0 o on o no o 


RaOImnC— «InFI/T*'OuTP<iTi PUNCH. T^Ptl0lTAPE60■INPUT^ 


RAOIANCE PROGRAtr 

UPWAH.I OH UOWNMAOD RAUIANT P0«ER 
blaCn HoOY CaLIR'JATION 

: - yggj Smith ^N0/0R OLD tables 

XATCR VAPOR aDOEh i>EC. 26. 1973 
PROFIi t« AOOEO JANUARY. 197» 


dOMMOn « (200. ^) . ^*r^ (200* ?> . T (200. 2). TO (200. 2> 

ComhOm uElUC(206.2) . *ODElu(2‘0.2)'. OELU^Z (200.21 
COmmOm XXOELU(200). XU£luC(2oO). XDELOZ(200» 

COmmOm sect (3)* All I (3). «/ C (3). THEtA (3) 

COMMO.r INELAG (307. Names (3'l) i INPur~ nO) 

COmmO.. TTh (2>.'A(KW (2». XPP (2). XTT (?) 

COMMOm IFmT (9) _ ; 

DImENcION IMU (^. 4) 

(llMtNcION WSS (200) 

C0MMO.i/RTW/XP( 200) .XI (200) .X y(200) ,TA(2 oO) 

COmmO. /head /ISTa. RuaTe, NO^TE. wlI* TRI-2. IFl. IF2. FNUl. FNU2. 

I'NHAGp. ONU. nPIm. IbRlfE 

COmmOm /XYZW /WX. wY. -2. Cw (9. 1?0) 

(iCMMO.. /XYZ /Ul). UY. OZ. CCO (9. (.3) 

CUMMO.. /AlFLIR /al (22V), Fltm (229) .' DI)HMY2 ' 

COmmOm /AuCaCw/' ALC(22V) .ALR(229) .AlOZ(»29) .0ZMXR(19) .0ZPHS(19) 
COmmo.i /TBANS /MAI. I*|aO. NtWA. TT (50) . TCOI. TCOO. NTCO. T TC (7 

14) ,Tj7I*ToZOJnTOZ.TTOZU1) " 

COmmOm /XWT /X (3). RT (3). XXI (6), XK? (6). 8B00Y (75) 

COmmOm /RaRAMT /M. NANGLES. Ul* NU?. Nil3» HK. LUQO. 

IPI ,LUf^.NoZSW ‘ 

COmmOm /waTERP /KHObS. FZW. .421* VLAH0 

COmmO». /OaTT /nCa. nCo. nwa. NWb* IPP. tSPRAU. MINOT. nPAER. PlOw 

1. PHIcH. CKNHG.’AL0.'ZtK0," CC22 *' Pup. POO. NKOwS. 'LCON^'tEVEtS'.- I r 

2**AIT. NbbRAU. LtVUP. L£Voo. TI. TO. TF . NLINES 


GImENcIUN NNN(6) 

EOUIV.LtNcE (fwMl. TWII). (T'WAO* TwiO). (NT**A. NTwl) 


c 


Data (»<Niveo. 0* o."i .“1 *“i) 

UATA ,*.7 * .bbSbSSbSSO* ,88b‘ibB88«9, .5955555556) 

UAjA ,X * - .774b9bb692. -.0. .774596669?) 

UATA (UEGTRAO .017453293* . TPI * 3.14)59265) 

UATA (ZERO = 0.) " _ 

UATA ‘(P « 3.74l3t-l2). (Q « r.4389) 

Data ,xxi ■ .ibJ . 12 . . 09 . .-.s, ,o7. . 07 ). (xk2 ■ 25 .. 19.. lo,. 1 


10., l7.. 27.) 

UATA(tHU» 8H WAitH ,bMYAP0R,7M CARBON. 7H010XI0E.8H TIIH. 
r3HUOk..?H fIl.jHTER) 

UATA (EPS « S.E-b) 

-data (A « 357.9110836). (DO ■ - 1 .226463094E*4) . TDl • - ^.4t02l5S 
U2E*1>. (02 i 2;i7l27bV6E-l) . (03 ■ - i .317730402E-4) 

IOaTE ( 1 ) «8H0U/MM/ YY 


I 
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nv the 


call cYbTEMC«30»NNN> 

NC A ■ t)oO 
NC0 B Hlu 

NWA B 6«lfl 

NW(:t B i^oo 

HSG B <'.U7E*t> / 9bu« 

— NOATf B IDATt (AUUMMYI 

l500 format <10Fb.2) 

NT«A r MThA - 2 

NTCO rNTCO - 2 

names (1) B ttHFMtUuENC 
names (2) B toMANliLtS 

names (3) B 6HFILTE« 

names (^) B JHC02 
names «S) b bHU>*BAKU 

names *(0> B HHOOWNbAHO 

names U) b /hPHOCeSS 
names (»> B OHSIUP 

names CaLIB • 

names (iO) B mhki*< oaia 

NAMES (11) B 7MSTrtTI0N' 

names <i2) » BMSPECTRAL 

names (i3) b tiHlH BElGri 

names (U) b bHOHlOMfNt 

names US) » 7MAEROSOL* 

Names (lb) B 7Hbl6NEL(. 
names (17) B bMlMbLES 

names (i8)‘ b «hnet temp 

NAMES (19) B 5MU20nE' 
names (20) * BMbATER VA 

NAmES (21 ) ”8~t>HPR0FItER 

IbRITr « 2 
ISPHAr B 2 

jpp ^ j 

NCOS- B 1 
NOZSw.l 

IFI B IF2 B IM " 

NUl B 0 
NU2 B 2280 

NU3 B 10 

ALO B 0 
AHI b 9u. 

LEVUP * LEVOO B 1 

PUP a lul2.8 

PUO * .1 

■ ■ WLl » 1. / 22b0* 

l'"L2 e ttMlNFlNlTY 
IIoINDt b 1 

NPaEP b 2 

NUNAMrS B 21 
NCAROc « 0 

IF (N^NAMES .LTi 1) (iO'TO“Ab.< 

UO lOx 1 B li NUNAmES 
loo iNFLAr. (I) B 0 

46S CONTImUE “ 

C****B**B**«*BB****B********.***************************** 

c read control CAROS 
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— c**** 

NX > >.U''<IT0H (NiiNMMt^* NAMtS, 1NPUT« 4Q) 

NCAHUc » hCAHllS ♦ 1 
WWITfi; UOt 13;»4llNHUl 

GO 10 t20ti» 215* 2j0f 2<»0* 2 5» 255t 3lSt 265f 215« 
1 165, Wo, 105, W5» 105, l<* , 150, 15S. 160, 305, 

2B0, 300* 
265), NX 

.1 

310, j 


; 

*200 

INFLAr (11 * 1 

UtCOUc.- (H, I5l0» Iw**!!! (31)N.l 






UtCODr (A, 1510* iNHUi 14))N ? 
UtCOOr (H, 1510* iNt'Of (SllNiii 


t . 

4- 



isio 

(•OHMAt (In) 





IF (N..2 .fcQ. 0) MU2 * 22B0 
IF (N.i3 ,E0. 0) MUJ * Iw 

f’nui - •'*Ul 






FN112 W'j? 




c**»* 

UNil s NUJ 





20b 

WLl « !•£•♦/ FN(i2 
IF (F,.Ol .EO. 0.)205* 210 
NL2 r aHlNFlNlTv 




i 

210 

GO TO W5 

»«L2 * l*tA / FNUl 

tNCOOi (8, 1512* I»iL2)*iL2 





\^U 

FOOMAt (F8.2) 
GO TO 1*5 



1 


— 215 

INFLAr (2) « 1 

OtCOUr (b, 1514* InPU( (2))A|0 
OtCOOr (H, 1514* iNPUf (3»)A<1 



1 

■ 1514 

FOPMAt (F8,C) 

If (a,iI .iu. 0.) GO 10 22<; 

IF (Awl .to.' AL(i) GO TO 225 





Cl a (AMI - ALO) • .5 
C2 a (AMI * ALU) 4 .5 

CC22 - PI • (SI*'* (oEGImmU • (AHl - ALO))) • • 2 



i 


Do 22' K a 1* 3 

IMFTA (6) a Cl • X (^) ♦ C2 

■? I OrGlHAO * (IHETA (t^) - A| 0) 




22u 

COST ; COS (T) 

4SC (k) a WT (K) • SIN (Y) * COSY • Cl * *10966 
SECT (K) a 1. /C05 (UEGTHAO • THETa (K)) 



1 

225 

NAWGLe-S ■ 3 
GO TO 145 
NANGLrS a 1 





aSC (i) a .invbb 
(META ( 1 ) a ALO 

5ECT (l> a 1. / COS (ObGTRAO • ALO) 




Cmmmt 

AlII (2) a AlII (J) -0. * • 

GO TO 145 




23y 

INFLAr. (3) a 1 
IFl a IixOuT (2) 
if 2 a INPI'T (J) 





call nUMMYS 

If (1».PUT (4) .to. JMNt»)?35, 145 




235 

HEAO ?5w0, FLTR 

. - « 

- 
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i 

I 

- -• call *.lnhrnt (Flih, iMU ii« 4 )) • 

Go TO !<»•» 

240 INFLA« U) ■ 1 

IF (I».PUT (2) .tu. 7H1 nCLUU£i NCOSw > 1 
if (If POT (21 7 *iLACLU0E< NCOSm > 7 

60 TO l45 ■■ ‘ ' 

C . UPWAHm _ 

24b INFLAr. (b) > 1 

"IF (If.PUT ( 2 ) ~.t0r5»iat0lN) ItECOOE T«»“751A* iNPUr (3))PUP 

IF (I^PUT (5> ;tu; otiLtVELS) OECOOf ( 8 , IblO. INPUT ( 6 »)LEVuP 
i> (PmP ,E0. 0>2b0. Ia5 

I50~ INFLAr. (5) i ff ‘ . 

60 TO Hb 

C — . UOwNW RO — 

2Ss INFLAr. < 6 ) « I" - ’ ’ , 

IF (ImPuT (2) .tu, btttitUlN) '>EC00E ( 8 * l5U» INPUT (3»)POO 
IF < 1 ..PUT (5) ,tu; 6 HL£vclS) UCCOOE ( 8 , ISlOt INPUT (6MLCV00 

— IF (Pro .eQ.''oi 2 bor'iAb ■ 

2o0 INFLAr. (t>) i 0 

60 TO 1*5 

Process* 

3ls M m M 

, c************************** SET Up to PRINT INPUT PARAMETERS 

FlLlNte»J 
NPA 6 E B 1 

IF (1.FLA6 (V» ,LU. I) 60 T'> 330 __ 

I'F (I. Flag (2o) iEuT II 60 *T ) 335 

i'F (I. FLAG ( 21 ) *LU^ I)' 60 tn 360 

Call rAL*»AT (x1 »~Fa I» 2 ) 

(iO TO I>*b ■ 

£•••••««•••••*••••••••••••• dB CAi IBRATE 

330 5-ALL rALMAT (Xli Fxl* 2) 

- GO TO 1*5 * 

C************************** bATER vapor LOOP 

335 XI « ,01 ‘ 

OElx * .0001 

Call tAlwaT (xi» Fxi* i) 

Call talmat (xl * utLx* Fxii, 2 ) 

FPXl ; (FXll - l-Xl) / OtLX 

X2 « vl - FXI / FPXI 

3A0 01 ■ »2 - Xl 

IF (X5 .LE. .0000) 60 TO 350“ 

Call tALMaT (x 2* Fx2t i) 

If (AbS (FX2) - EPb) 33bt 355 . 3*5 

3*6 Call rALuAT (x2 * 0ELX,*Fx22. 2) 

PPX2 i (FX22 - FX2) / OtLX 

AI B <2. • FPX2"* FPXI - 3. • (FX2 - Fxf) / 01) / 01 

- - UI B rX2 / FPx2 

X3 B y2 - ul * 11. ♦ Ul • AI / FPX2) 

Xl B »2 " _ _ _ 

Pjjj ^ ^ 

X2 B »3 

FPXI « FPX2 _ 

GO TO 3-*0 

350 PRINT ls30 

1530 format (//*v IS LESS Than or equal TO .O OOO* //) 
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tALL 

call 

Call 

l«ll 

Call 

Call 

00 TO 


rALWAT 
rALwAT 
rALWAT 
rALWAT 
rALMAT 
rALKAT 
IA 5 ~ 


(.uuit >-Aif \) 
(.OUlfOf t>Aif 1 ) 

(.UUdt FXl. n 

(. 0 V 7 O* FAl.'l) 
(. 02 Sf fXIi 1 ) 
(. 090 t Fxit n 


3Ss T ■ T A ( 1 ) 

Prt » *H ( 1 ) 

ESLN ; A ♦ 00 / T ♦ or • AtOTi (Trv“ Toy»”T-« 021 ■• T 

LS « ^ 0 . • EXh (£Sln) 

WS ■ • tS / - ES> 

RH B Cz / US 

PRINT IS32. HH 

Format <///• RtLAriVt HUMIOtTY 
60 


1532 

C***< 


-•/» 


3o0 


260 

1516 


1518 

270 

A 7 o 

1520 


1522 
C 


•.F7,?t* 

TO 145 

>«•••••••••**•«••••••* temperature profiler 

Call n«OFILE(M»NPTwiALOlNU3) 

GO TO 145 ~~ 

_ STOP !-,H EOF 

HE* I Nr 1 0 

PRINT IS16 

format (-ilist of input conTvOl cards F0LL0*»<S*/) 
IF (NrAROS .CT. 1) 60 TO a7o 

00 27= 1b I ♦ -NCARoS' 

read (10. 1S2A) Input 
PH’NT i>»iB. Input 

Format Tix.loAB)-^^ 

CONTI.tUt 

CONTIhUL 

format <1H) 

if (ImPUT (I) .to. TrIPHOCESS) PRINT 152^ 

PRINT 1P22 

format (*0ENU OF HUN*» 

call fXIT 


275 


OB CA, IBrATE 

INFLAr. (y) * 1 - 

INFLAft <19) » INFLaO (20) b 
uEcoOc (8. 1S14I Input (3)) n 

UtCOUr" (8. IS1 <»(-InPoT-<a) ) TF" 

OtCUUr (Ht I9I4* InPoT ( 5 )) To 
60 TO l**b 


P T M OATA- 


280 


290 

1524 


iNFLAr. do) ■ I 

ENCODf (72. 1524* IFMT) (INPUT (1 ) *I«2« Iq) 

M s 0 

HEAD Ti924, INPUi 
format (10A8) 

NCAHDe a NCARDS ♦ 1 

WRITE do* 1624)INPUT 

IF (I. POT d) .to. 7RPIW END) GO TO 295 
M . M ♦ 1 ' * — - — 


ENCOOf ( 80 . 152 a. INPUT) input 
UECOOf ( 80 . iFMd lNPuf)XP 
TA (H) * XT (M)**’ 273 .i 6 
T B T A (H) 

PR B yP (M) 


(m). XT (M)t am (M) 
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ESLN , A ♦ 00 / T * Ui • ALOr, CT) ♦ (03 • T • 02» • T 

bS a fO> • bAP (tSuNi 

aSS («•» ■ 622 . • bS / (PR - F!§) _ _ 

IF (WcS (M) .LT« 0 .» (M) a 0 . 


(»0 TO 2vO 
NPTP • M 


1526 

Print i» 26 t (*p (I)*‘x« (d« xw (d, ass <i)t i ■ u nptp> 

Format (lHlt 5 X»lMP, 9 X*ItlT. 9 X.lHa» 7 X,aP ^ATa/( 2 F 10 . 2 . 2 F 10 . 6 » ) 


00 TO 1 A 5 


C — 

300 

iNFLAft ( 11 ) a 

1 


ISTA ; input 

MOATE ■ INPUT 
00 TO 1 A 5 

( 2 ) 

( 3 ) 

310 

ALO a 0 * 
ISPWAr, a 1 



00 TO 225 ■ 


169 

IPP a *2 


INFLAa ( 13 ) ■* 
GO TO l <*5 


^l 7 o 

175 

16 U 

INFLAr. ( 14 ) a 
IF (I.PUT (A) 
IbRITr ■ 1 

1 - - 1 

.tu. 7 H 1 NCluOE; 175 . 180 

GO TO 14 S ■ 
IbRITr « 2 
00 TO 145 


-- C 

185 

3 

1506 

iNFLAr. ( 15 ) * 
NPAER a 1 
tNCOUr ( 32 , 3.1 
t URMAt (<»AH) 
UbCOUc ( 32 . 1 
format (EB* 2 . 
IF (AoSA .EQ. 
CKNKG a ABSK 

At R 05 U ••••••••. rnmmmmm ••••••••••-»—•• • 

1 

NPUT( 2 ) ) (INPuKT) .Ia 2 , 5 ) 

506 . InPUI ( 2 ))a 6 SK. PLOW. PMI 6 H 

3 F 0 . 1 ) 

O.'I 00 TO 190 
• kSG 

- - ... - 

GO TO 1^5 


190 

INFLAr. ( 15 ) ■ 

0 


N«-'AER » 2 



GO TO 1*»5 


19-3 

INFLAr ( 16 ) ■ 

1 _ 


IWINOt a 2 


^ 

GO TO 1*»5 



IF (WPJT ( 2 ) .Nt. ‘»«L16T) Gi TO llO 
<JalL wL'^PRNT (ALt'lHU (1. D) _ 

Call wLnprnt (alc* imo (i. 2 >) 

Call wLNPRnT (Alw* IHO ( 1 . 3'^) 

Call . lnphnT (FCth, Iiiu a, 4 )] 

— Call twrteR (a» b) 

Call tCmkhon (Af b> 

60 to U5 
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UO~AP (U**UT (2) .tUr «H*»ATER Va)ll5t 120 
11a HEAO T5J2. na, Nb" 

1502 format (2I5) 

HEAD TbJA* <AL~<NN)»-NN m'NA, N8) 

150<* format (IOFB.2) 

Call wLi'*RpNT (AL» IMU (if D) 

GO TO 1<*5 

12o AF (1..R0T (2> .tJ. brtCAHBON .i)l25t 130 
125 head i502f NA, Nb' 

READ TSOAi (ALC~(NN) I ■ N I NB) 

Call wLnrrnt (ALCt imo (it 2>) 

GO TO US 

-l'3H ■ iF df.PJT' i JSf-30S 

135 «EaO i5c.2f NAf'Nb' 

HEAD ibOAf (ALW (NN) I NN s N^f NB) 

Call wLnprnt (alk# i«ir'-tif-3r>r 

GO TO US 

C C _ NET TfMPERATuRE 

140' iNFLAr. (lb)~»’l 
GO TO 14S 

C . OZONE 

ISO " INFLArt (19) » 1 - ” ~ ^ 

IF ( 1NdU((2> ^nInCLUUI- ) NO/Swsl 

IF (INaUr (2) itO.fMLACLUyt) NO .Sw = 2 

(io TO U5 " 

c- water vapor — 

155 iNFLAr, (19) ■ INFLAG (9) ■ o 

- - INFLAr. (20) » 1 

EncoOf(2a, I t Input (3) ) Input (3) t input ( 4) t tnPUT (S) 
1 fORMAT(JAtt) 

UECODr (24t ISO? t IN^T— (3) ) vLAMBr RNOBS 

1506 Format (FBtOtFiP.o) 

GO TO US 

c — - profiler — 

Ibo INFLAr, (2D ■ 1 

iNFLArt (9) ■ INFLA(i (19) » InFlAO (gO) « 0 
GO TO US " 

c— . — Incorrect input control 

305 PRINT i52bt Input 

1528- format (•OlNCORHECT-INPUT-COnTROL«/lH<TtrOAal 

60 TO U5 


0 0,00 *-00000 lo o oo o o oooooooooooooo 


I 


fUNCTtON MONITOH <NX* NAMES# INPUT# t.Ul 

////////////////////////////'^-'^^ f ^ 

TrilS Function heads one CAHU FHOm logical UmIT lu. the first eight 
COLUMNS crNiAiN A sihing'o^ a hollemHh chaoactehs. This string of 
CHARACTEPa is COMPAHtO wIlH The i IST of R-CwAHaCTER STHINOS stored TN-- 
the array names. if a maIch is found, the value of monitor is Set 
EOUAL TO the SUaSCRIPT OF The matching name in The list. a maximum 
OF NX NAMrS will BE 'SEARCMtO; — Tf NQ MaTCHTNQ NAME IS FOUND# ThE 
VmLUE of mO<^ITOR kill dE^bEl TO X ♦ I , IF ThE CARD WHICH WAS 

READ* IS Am END-OF-FlLE CARoJ THE VALUE OF MONITOR IS SET TO NX ♦ 2 . 

IN the aRpAY input WILl^8E■^^ETUR^ ED ten WOroS OF HOLLERITH' CHARACTERS 
ReAO'FROM COLUMNS l-bO OF'IhE CAr(), ThE ARRAY InPUT MaT 8E USeO WITh 
A DECIDE STA TEM ENT. 

oavio l# obitts 

JUNE 1968 


//////////////// / ///////////^//-'^/^ 

D iMTWmJ^NAMES m r INPUT TTU) ~ 


f 0 i^'M 5 ~T s 


FORMAf“TroA8I 


. - _ - - read one CaRO with a control BOrO in columns 1 -S. 

head (LU# ISUOIINPUT* 

■ If (EoF( lu) ) 100# i05 
c--“---AN end-Of-file Card as read. 
loo monitor « NX * ‘ ■ 


RcTUHm ■ ■ 

los name , INPUT (I) 

IF «Nv .LT. I) Go TO 120 


110 

00 lln I « 1# NX 

IF (N/>ME .EQ. NAMES (I))115# 110 

COnTImUE 


120 

C - - 

CONTImUE - 

MOnIToR * NX't l' 

• # 

C - • 
115 

■ RETURm 
MOnIToR ■ 1 ‘ 


- KETURm 

End 


REPHODUCIBILITY OF THE 
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— — SUHKOitTlNf CmlwmT (wu* "ODE) 

COMMOk' •» (200» 2)* pH (200* 2)» T (200* 2) » TO (200» 2) 

COmmUki UELUC <200 *2) * OOtLtl(2'0*2) • DEtljnZ (200 *2) 

COMMO*. AXOELO(2VO). XOtLUC(2»0). XOELOZ(200> 

CQmhO*' sect (3)* AIIl ^3)* (3)* TMETA (3) 

CUMMOm INFLAW (JO)* names (3a>* INPUT MO) 

COmmO»> ITh ( 2 )*'x*«iii( (2)* XPP (2)* XTT (?) 

COmmOm IFMT (9) 

DImFNcION FNET ( 200 * 2) •"VNEr (200) 


C0MMnu/MTW/XP(2U0) *XT (2uO) ,X»<(200) *TA(2rtO) 
• -yLjr“Ar\ t r<T K . An Air. M,.ATr. Uii."tUL 


COmMOu '/head /I^TA. AdAIE. N.'ATE* WlI* 
INPAGri DNU* NPIh* IdHlTE 
COmmOm /XYZw /wX* wY* CH (9, 120) _ 

COMMOm /XYZ /ox* OY* OZ* CCD <9* 63) 

COmmOm /AlFLIP /AL (229), FLT« (229), niiMMY2 


TWL 2 » TR»- ir f* FW UI, f w ur ,— 


COmmOm /AiFLIP /AL (229), FLTR (229), niiMMT2 

COmmOm /AlCaC^/' ALC(229) ,AL*m229) ,AlOZ(?29) ,0ZMXP(19) ,0ZPHS(19) 
COMMOm /TRANS /[tXAI, T*(AO, NTWA, TT (50)*”TCOI, TCOO* NTCO, TTC »7 


XK? (6) * BBOOY 

rNil3*~MK,‘LU00, 


lA) ,T07 I*TPZI)*nTUZ,TT0Z(‘*1) 

COmmOm /xwT /X (3), '*T I3),.x*(l (6) f 
COmHO,' /PAPAMT /Mr NANGLES»“ Ul,'NU2 
IPI ,LUn,No2SW 

COMMOm /waTEPP /HNObS, FZW, -2T , VlaM0 

CummO*. /UATT /uCa, nCcJ, 'nWA,“NW6,-IPP,~tSPRAD, IWINOT* 
N pMIftM, CKNP6,~AL0 ,“'ZEm 0, CC22, PUp, pOO, NKORS, LCON, 

2«AIT, NoHRAO, LtVUP, l£'^00* TI* to, TF , NLINES 

- COMMO*' /value /All ^ ' 

Dimension nnn(6) 


( 75 ) 

NCOSW, "Or 


npaer* plow 
levees, IR 


C 

c 

c 


-•£QUlVALENCE-(TWAr,— Tirt,iri-Tr^«l“TWtDT,-(NTWA,^ NTWI) 
DATA (mNN« 0,0*0 *-l *-l »-P. 


lOU 

730 

10a 


noza; *970 

NOZBl.lvJBO _ 

NOZA2b1o8o 
NOZhZsUjO 
CALL rYSTEMC(30»NNf^) 
GAmmA a ,b5 
IF (ImFLAG (20) *£0, 
IF (M ,LT. D GO'TO* 
bo lOr 1 a 1* 

XW (1) a wO * 
COnTImUE 
CONTImUE 
Xinw (?) < 


0) 

730 


GO TO 105 


/ XP (1)) • 


VLAmB 


XmW 

XPP 

XPP 

XTT 

XTT 

TTB 

Ttb 


( 7 ) 
(?) 
( 7 ) 
(i ) 
( 7 ) 
(?) 
( 7 ) 


XW 

XW 

XP 

XP 

XT 

XT 

TA 

TA 


( 1 ) 

(M) 

(1) 

(M) 

( 1 ) 

(M) 

(1) 

(M> 


0 £r'lN£ u Values 


100 


1 


I 

1 

1 

J 


J 

] 

] 

J 

J 


MM S M • 1 ' ■■ ! ! 

iF *mm .lT. n teo ro 7Jb •• j j 

00 11a 1 ■ l» mm = 

Nt1*1 

AtJSMM ■ AUS (XP (I) - *p <NJ) i , 

AVPWS. (XP ( 1 ) *XP In) ) ••ti 

WMAH , .5 • m ♦ Aw (NM i 

iF(AVc-Hi.LT.OZpHSU) ) 00 TO 

b2ttAH;oZMx«(i) _ f 

GO TO 17 ■ 

16 lF(AVrHls.6T.0ZPMS(l0) ) GO TO 18 ' i 

0Z«AHe0ZMXH(19) ' ‘ 1 

GO TO 1 7 I i 

18 00 IS J0«1 ,1V ; 

1F(AVdHS.lT. 0ZPMS(J0)) 00 TO 15 i 

— *6zBAH;((AVPH5-ozpHs(jo))*<oZMAR(jo)'ioZMxRrjo'-ii)7/iozp«?rcJO-n- 

•OZPHS( JO) ) ♦OZmXH( JO- i) ' I I 

GO TO 1 7 ' J. I 

l5 " CuMTl*Ut - -- , 

17 XOtLO? ( 1 ) *ArtSPP*OZdAW*»00102 i 

AOELO? < 1 ) *AOSPP*OZbAK*rw0102 ‘.A ^ i 

AxoEUn tl) * -bah * AaSpp ^”0010204 ‘ | 1 

llO AOELUr <I) » AdSPP * .2*»8 ' I i 

73S CONTlK.Ut - -- ^ *!_ I ^ 

c*#***i***»*****************v****************************************** ® i 

c HEVEHSE P«T,w» '.NO U ARRAYS i 

C*##**»*<^«»*#o«***********^********«**^******^**^****^**^****** ” * ^ } 

NPtmi ■ NPTG _ j; ^ 

NPTPl « NPTw ♦ 1 t 1 

OOELU (NPTW* 1) « OUtLU (If 2) » 0. 

- IF (NoTw .LT. 1) 60 to 740 - - - *« . . 

bo 12' I » 1* NPlw 

jj « uPTPl - I ^ 

W (I) » w (JJ ♦ 2uO) ■ Ak (J«) - . - i ; 

PH (1) * P« (JJ ♦ 200) » XP (JJ) i 

T (I) » T (JJ ♦ 200) « lA (J )) . 

to (1> « TO (JJ ♦ 200) ■ XT (JJ)' ' ? 

If (I .to, 1 ) GO TO iiG 

t (JJ ♦ 200) a I d - 1) ■ »S • (T (I - 1) ♦ T (I) ) 

If (I .to. NPTWl Go'10'l20 
IIS JJ « *iPfW - I 

OELUC (1) « OELOC (JJ ♦ 200)_* XOELUC (jJ) 

~ OOELU (1) « OUELO (jj ♦ 200) * AXOELU (jj> I 

0CLUU7 ( I > *0£(-UQ4 ( JJ*200) sxOEi OZ ( JJ) | 

120 COnTI'.UE _ I 

’7a0 COnTImUE •• 

12s If (IvPOT (2) .to. SMSMITh) l->9» us 

130 ASSIGm 305 TO N§MITH^^ j 

60 TO l-*0 ' " ' * I ; 

13s ASSIGN 310 TO NGMITH i I 

.'oSlGf 355 TO IPIM _ J 

assign A05 TO iPTu J 

140 IF (I*. FLAG (V) .£0. I)lW5. mo ^ 

14S <»0 TO S7o ' ‘ 1 

15o "IF (U.FLAG (5) .EQ, 1) 155* l70' i ^ 

15s LUA » 1 ? 

IF (U.FLAG (ft) » £Q* dl<*0» j 

/ " \ I 
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"lou* LUB » 2 

60 TO 17b 

16b LU8 « I 

— GO TO W5 - - 

170 LUA B 

LUB « 2 

— 17b IF (L.iB ,lT. LUA) 60 TO 745“ 

liO bS«? LUO “ LUA* LUB 

LUOO i 200 * (LOU - 1) 

MKI B 1 • " 

60 TO (180* 18b) • LUO 

180 LtVELc « lEVUP ■ 

MKF B NbUBS (PUPr PR“tlTl) »~WPT»()— ' 

60 TO ivo 

18b LEVELS » LEVOO 

MKF B NbUBS (POO* PH (1* ?) «-NPTW) 

190 MKF B MKF 

Mkl b MKF - 1 

IF (li.FLAG (21)“.EQn) 60 T.V 220 

IF ((uPMGE .to. i) ;anO. (NLIMES .EO. 0))195» 215 
19b if (IrFLAG (2o) tEO^ U> 60 TO 200 

- l'aul mEaowaT (LUO) 

PHINT IbOO* VLAMd 

1500 format (90X**LAMHOA ■ •*F6.2) 

ULINEc B nlINES ♦'! 

60 TO 21b 

200 00 TO (210* 20b)* IPP 

205 PHINT lb02* ALO*'NOATE * 

1502 format (1h1*VX**ANGl£ ■ •»F5.1*40X,bCOmpuTEP **A8) 

nlINEc * 1 ' * 

GO TO 215 

2ly Call uEmDING (LUU) 

NLINEc ■ NLINES ♦ J. _ _ 

2lb IF (1».FLAG (20) .EO." 0) GO TfV225' - “ - - 

220 MKN b 1 

mkl b mm 

60 TO 285 

22b 60 TO (230* 245)* IbHITE 

230 ASS16*i 235 TO IbLAQK 

• - temp , - 770. 

00 2**a J b 1* 7b 

60 TO 605 

235 BHOOY (J) B buMO — 

240 temp c TEmP *5. ^ ^ 

levelsloop -- 


24b ML a T 

MkO * 0 
MKN B MKF - 1 

250 temp i fT8 (LuO) - 2 / 3.16 
ASSIGk- 255 TO IttLACH 
GO TO 685 

255 60 TO (260* ?65)* IBHUt 

■ 260 AII-b SUM6 ““ “ 

UCX a 0. 

UUX B 0. 



oepkoducibilotofthe 
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ASSIGoi TO NbKITt 

GO TO bOn , 

If (1»FLA6 (iO) ttU. u) 60 T) 270 

MK s mKO .1 ' - • - 

FNET jMk. LUOJ ■ SUH6 • PI * l,E*A / 69P* 

if (.UFLA6 (i3> *^0. 1) GO 260 

PrtINT IblBt 2EROt SUrt6»-ZER0. ZERO' 

SUM? T GUMo * C022 

print U20* App (LUU)t ATT (| UO) ♦ XWW (lUO) t ZEROt S0M7 
IF (MSL' ,LT. MKIT-G0^0-75(r ' 

MKsMKt 

Mkn ■ MKF - MK 

^ angular loop 

•«••••••••«•••••*••••••••••••••••••••••••••*••***< 

K»l ■ V 

CTT « SECT (K) 

MKI b mnn ' 

MM s f.'KL 


frequency loop 


SUM2 r 60^3 * SUMS * SUM? * a; 

SUMTA b SUMTe B 0, _ 

IF (N.i2 ,lT.‘ nuU’gO to 760 
UO b6e 1 8 NUl* NU^it NU3 

fact ; 1. _ _ __ _ 

If ((t •£0» nuD »or» ct .eq. 

UNU S I 

IF (I»'FLaG <I3) .EU. 0) GO Tm 


NU2J) FACT 


IF (I»FLaG <I3) .EU. Tm^295 . 

print ISOAt UNU * " ■■ 

format (//•OPHESSUHE*»I'JA,*TrMP#«l2x»*OT/D(LN(pn*.lOXt«FREOOENCY 
l« *.Fa. J ) ' _ 

Punch i306* isTAt muate# unu 

format (13XtilA8*«F*9R£QuENCY.il ■ * ,F 4 . 0 . * < l/*9CM4-l ) •) 

MKI » MAI . ' 


UNU3 c UNU • * J 
CX s ,.X s UUA B UCX 
AUX s XCX * 0. 

Pi* s tR * PC B IC'b 

O2x*XpZasuO2Xb0» 

P02bTpZ*0, 

SUmTA * 0. 

ipl » 1. 

Go TO NSMITH* (UOSt 
ASSIGfi J55 TO IPTM 
ASSIGm A05 TO IPTO 
IF ((iiNU ,GE. 2u2.b) 
IF ((.iNU ,GE. 507. S) 
fTpLUr • TO (MM'*"ii 
SUMS r Sum b u. 


b) .AND. (iiNU .LE. 
5) Iano. (mNU .LE. 
1. LUU) ♦ :;73.16 


797.5) > ASSIGN 390 TO iPTh 

817.5) ) ASSIGN A40 TO iPTO 


height loop 


IF (Mw .LT. MKD go 10 765 
00 55r J B MKI. MM 

jcluo ■ j ♦ Luoy 


S 



II 

Au* ■ AMx ♦ untui <J* »-oi») 

XCx a xcx » UfcLUC (j* LOU) 

XOZXa.O^X«l)LLuOA < J*LOU) | 

••••*•••• 

c CAUCULArt tFFtCTiVE PRESSuiJE AWO TEMPERATURE ' 

* 

OUCX . uElUC (J* LUU) • CTT * I , 

UUUX . OUfl.U Ut LUU) • CTT | 

UUOZXrOtLUOZ ( J*LUU) «CT f • ' 

Pri a .b • (PH (j* LUOr* PR VJ ♦ 1‘, LUD)) 

31s UUX a UUX * UUUX 

PW a nw * PH « UUUX ^ I 

— TW a T* ♦ T (J. LUU) * UUUX ■ 1 

3^«) a Pa / UUX 

32b tEa a Ta / Uux 

330 UCX a UCX ♦ UiJCX 

PC a nC ♦ PH • OUCX 

tc a tC ♦ T (J» LUU) * OUCX 

‘33b PEC a PC / UCX - • 

340 tEC a TO / ULX 
00ZXaMOZX«UUUZX 

POZaPnZ*PH**OU()ZX 

tOZ*ToZ*T < J»LoO) •UUOZ a 
PEOZ aoOZ/uOZX 

TCOZaTOZ/liOZx 

PRftT , PH * .OOi 

UX a ,,X * UUUX • (PHAN • • -AMMA ^ 

CX a r* ♦ OUCX ^ (pHAf) • • - - '• 

0ZxaO-rX*UUOZx* (PRAT) **.Hb 
J4b IE <U»)jH«.* J»5* JSO ■ ' 

350 Go To IPTH» «35»» 3<#0) -■ ' -- 

355 (hat b 293. / ( (J» LUO) 

If <trcT.le.o.) uo to jso 

A1 a cLOGlO (TRAT) 

35b A2*TR»T-). 



C - LOG(L) FOR water VAPUH • ~ - 

* 

IF (Uy .LE. 0.) 'GO TO Jb5 

IF (I ,bT, 95o) 00 TO 3b0 ' 

UOIF T UNU 
ACoEF a ,<»8E-5 

(»0 TO 3T5 ■ ■ “■ ■ 

360 OUIF r UNU - )5Vb. 

IF <UnIF)365* J/Ot HO 

365 ACOEF a 2.AE-5 ' ' 

GO TO 3T5 

37o ACOEF a 1.7SE-5 
375 WW a T • .1 

IW a k.R 

ALIll a AL (la ♦ 1) ♦ (aW - f ) * (AL (Ta ♦ 2) - AL (la ♦ 1)) 

AAlOG « ALIII ♦ Ai - ACULF • A2 • UOIF • • 2 


C TAU F M aATER VaoOR 



ULOGa. a ALOblO <UX) 

ALU a ULOGWA * AALOG 


104 

!1 


i 



MW > (ALU - / FmAO ■■ 

IM B irfW 

If (lu .LT, 0) GO TO 3bb 


i 

! 


IF (lu .6T. NTMA) SO TO 3g0 

t2 ■ tT'UW * 1) ♦ (4m - JM> • (TT (IW * 2) - TT 
If (T? .ST. 1.) ’t2 ■ 1. 

(IW 

* 1)» , 

3B0 

IF (T-» .ST. l.E-10> SO (0 39g — 
TP2 b 0. 

SO TO 400 



3«5“ 

T2 » T. 

00 TO J9b 



~39o~ 

IF (UitX .lE. 0.) SO TO 385 ' 

Ux B .LUS (Ouxl ' 

MY B aLOS (Htw • .OUl) 


i 


MZ b aLOS (TLm V ifJtlty) 
T2 b tAUM (UIXU)‘ 


9 

39b 

40U 

TP2 B T2 

SO TO (400t 4301 ♦ nCOSm 
SO TO IPTD. (40i>» 440) 



40b 

C 

IF (Cy .UE. 0.) SO TO 435 ^ 

IF (I .LT. NCA) so to 4Jb 
IF (I .ST. NL8) so to 43b 

L06(L> FOR Carbon dioxide 

«#• 


410 

41b 

42u 

UOIF « UNO - 66/. 

IF (i)nIF)4lOt 41b* 4ib 
ACOEF ■ 4.6E-4 
GO TO 420 
ACoEF ■ 3.4E-4 

MM s j • . 1 


- 

c 

I M B vW 

ALClIt » ALC <Im ♦ 1) ♦ (MW - iM) • <Alc (H ♦ 2) 
AClOS b ALCIII ♦ A1 - ACOEF .A? * uOlF * • 2 

Tau F(iM CARBON OtOXIOE 


alC (IW ♦ i)) 

ULOGCo » ALOGlO (CX) 

AClU s UL06CO ♦ ACUOG 
WW X ,aCLU - TCUI) /’TCUD 


I M B wM 

IF (Iw .LT. 0) SO TO 4ib 
IF (Iw .ST. NTCO) SO TO 425 




T2C e TTC (Im ♦ 1) ♦ (mm -"1^1 4 (TTC (T** - 

IF (T^C .GT. 1. ) T2C * 1. 

IF (TpC .ST. lit-lo) GO TO 4<.5' 

TTC 

(Itf ♦ IT) [ 

42b 

430 

IP2 » 0. 
Go TO 4C40 
UCX s U. 


t 

43b 

T2C » 1. 
Go TO 4*.5 


-y 

•t 





44o 

IF (OrX .LE. 0.) GO 10 435 
(ix X aLOS (OCX) ' 

OY X aLOG (PtC • .001) 


■* 

1 


■« 


DRIGINAIJ PAGE IS 
OF POOR QUALITY 
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Ul w «LUO (TtC / 

T^C ■ tmuco <unO) 

C _ 

♦ ^P^ • TH 2 • T 2 C 

Go TO (1 * 2 ) tNOZ^w 
I lF(OZy.LE.O.» 60 To J 

T.tJOZAl) 00 TO 3 

IE (I.rT.NOZU<^) 00 TO 3 
IFd.l T.NOZBl) 00 TO •* 

- • IFd.cE.NOZA^r OO TOTST 

60 TO 3 


— t.06(LV F0« OZONC 


I* U0lE«i.NO-1043. • ' ■ 

ifniOTF) 6*7*7 — 

5 0UIF>mNO-1110. 

1F(U0»F) e*9*9 

— a acoef;i<*.e-a 

60 TO lu 
9 acoef;<*o,e-a 

00 TO 13 : — 

6 acoef;<*.e-a 
60 TO lU 

-~7 acoef;ia.oe-a 

iO 

ALOZlT*ALOZdw*l)*(»i*»-l«i)*(A| OZ(lw*2)-AL0Zd*<*l») " 

AOZLO«»»*LOZI 1 *Ai-AC0tF*A2*O0TF*«2 


— Tau fob ozone — 

#*#•••*«***•••#•*•*#••*******•*••«*•••••**••#•*••••*•••••••••* 

UL 0607 «Al 06 iO (OZX) 

AOZLU;ULOGOZ*AOZLOO - ' ' 

WW*(Ar.ZUU-TOZl)/TOZO 

Imsmw 

IFd«.LT, 0 > 60 10 3 : 

iFdw oi.ntoz) Go to il 

T20Z*TTOZdw*n * (Wn-U) •(TT07dW*2)-TT0Z<I«<*l) ) 

lF(T 2 r.Z. 6 T.l.) ‘tZOZsl. — - — 

IF (T 2 rZ. 6 T.l.E-l 0 ) 60 "f 0 12 

11 tP 2 «o^ 

Go To AUO ■ 

2 uozx>; . 

3 T 20 ZxT. 

12 iP 2 »Tn 2 «T 20 Z 

Go TO (<* 6 b* 450 )* IWlNUl 

c . ^-'oignEuI- window 

450 IF d - 700)475* 455* 435 

453 IF d - 1200)460* 47b*‘475 
460 w'w X ^01 4 (1 >*7UU) 

I* X uW ~ “ 

MJ X wM - IW 

VAlKI X XKl dw * 1 ) ♦ WO • (XKl (IW ♦ ?) - XKl (Jw ♦ 1 )) 
VALK 2 X XK 2 (IW ♦ 1 ) ♦'WD • »XK 2 (Iw * 3 )" 4 -XK 2 (IW*'!)! 
PdAR ; .5 4 (PR (J* LUU) ♦ >., (j ♦ 1, L'lO)) 

OtLP i abS (PR IJ ♦ r*'LUD) - PH (J, Lun)) 
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D 

D 


D 


fl 


i 


O O O I 1 o o o 


* 


TrtAH , I (Jt LUU> — 

WriAH r *SE-3 • (m (J * 1* LUrt) « <* (Jt |U0>) 

TmlN , oHAR • PdAR • UtLP • (VALKl • TtR (LUO) / TBaR * VAUK2 • «B 

lAR / ,622) / XPP (LUU)’/ 980. — 

SWRTH*» SUHTo ♦*FMlN 
T2m > Eap ( - SUHId) 

(i(> TO A«0 - - - _ . — 

C 6TAN0«R0 window 

Abd IF (1 .LT. N«A) go to 4/5 

- 5Q ^0 4/5 - * ■ 

t06(Lt FOR WINDOW 

If (Uy .LE. «.) (»0 TO 475 
ULOGW<t ■ ALOtolO (UX) 

— ww . r* M 

I w ■ biW 

ALwIIt • ALW (IW ♦ i) ♦ (ww . IW) • (Alw (IR ♦ 2> - AlW (IW ♦ 1)) 

AWLU > ULOGWA ♦ AUWlll 

WW ■ (AWLU - TWlI) / Twio 

Iw ■ b>W 

IF (Ii.. .LTt M> Gy TO A75 
IF (Iw .‘GT. NTwi) GO TU 470 

■ ‘ 'TAu F .R »(IND0 W ' 

T2 ik ■ TT (IW ♦ 1) * (<*w - IW) • (TT (IW ♦ 2» - TT (IW * D) 

If (T5w ,oT. l.)‘TiiW ■ X, 

IF (Tpt. .GT; l.fc“10) GO TO 4 m0 

♦ 7o 1P2 « 0 * ’ _ _ 

5o TO Ado ' " -- - - 

475 T2w b 1. 

480 IP2 « TP2 • T2w 

If (J - MKl)G)5t 5ibt 4d5 

485 GO TO (490t Gl5)t NPAEK 
490 IF (Po (J. LUO) PL(iw)i>15t fc95i 495 

49a IF (Pd (J, LOO) - pMlGM)5ov,t SOOt 515 

c — . aeros u calculation — —————— 

500 OtLP * AdS (PR (J ♦ it LUO) - PR (J, LUO) ) 

PUAH . .5 * (PR (J * It LiiO)“* PR (jrLnOn 
SomTa b SuMTa * T (Jt LUO) • OELP / PBaR 
iPAfcR B £xP ( - CknRG • SUMTft) 

If (Tpaeh - .oi)5osr bo5t 5i- 

505 MM s I - I ■ 

GO TO 290 

510 *P2 B TP2 • rPAtR 

515 OT B jPd - TPl 
TPl 8 TR2 

Go TO 025t 520) t IPP ■ 

520 PdAH s .5 • (PR (Jt LUU) ♦ P.i (J ♦ 1, LUO)) 

UTdLNo b pbAK • OT /{pH (Jt LUO) - PR (J ♦ It LUO)) 

PRINT l308t PBAHi'T (jJ LUU)‘. OTDLNP 
Punch loiot Pbah, utulnp 
150d Format <lM0tF8.2tFi4,2tF23,5) 

1510 format (F7.ltF9.4) 

GO TO 550 


— Q - — PLANCk FUNCTION 

•••••••••••••••••••••••••• 

52S r • UNu / T (JCtUO) 

WW • T • .1 • • • • - - 

FLTHItI ■ FLIH »!*• ♦ 1) ♦ <«<w - IW) • (FLTR (1M *2^- FLTR (1 m ♦ 

ID) ' ‘ ■ 

IF (PwHl ,LT. 6V0«» A1 ■ P • UNU3 / (Pi • <EXP (PwRD - l»D 

530' IF (J .LT. MMJ (io’TO SO" “ ~~ 

PWRZ ^ U • UNU / TTtfLUU 
IF (P..>H^ .G£. 6Vo;> (iu'lO 53 k 

AIO » P • UNU3 f IPI V ‘EXP (P«R2) - T.D 

(»0 TO 5A0 

53b *10 ■ 0» 

5au SUM3 ; SUm 3 ♦ AIO •“ 1P2 • FLrRIII •TACT ' 

54b SUM « SUM • AI • 01 « FLTRlIl • EACT 

55U CONTlnUt . . - 

7b5 continue ~ ' . *• 

C ■ * End of MtlGMT LOOP 



(»0 TO (355* bbS)* IPP 
55b SUM2 ; SUM2 ♦ SOM 
• ■ 00 TO (b60* 555) ♦"ISPRATT 

5b0 *11 » SUMS - SUM 

Mil « All / FaCI 

PHINT 1512. UNUJ Xn 

1512 format (2F20.8) 

56b COMTl>.Ut 

760 C0NTI*.Ut “ "■ / 

C ■ End OF FRtOUENCY LOOP 



(»0 TO 070* 040). IHP 
570 *11 * (5Um3 - SU.M2) * UNU 

— IF ((.Nt-LAG (20) .EU.' U) .AN' ,-(lNFLA6~72irTEO.~0n GO TO 5?r 

Fv»Z « All - HNOaS" 

If (IvFLAG (^1) *EU. I) 60 To 580 

— WP ■ rfO • 1000. ' * ■■ 

Go To (575* 360). MUUE 

57a PAINT l3l4. »P*'AII. UUX. F27 

1514 format (lH0.M5.2.t2*.5*E21.6*El6.4) ^ 

580 HtTUH*. 

58b IP (l.FLAG (18) *EU. O) 60 To 590 _ 

Ml > *<l • 1 

FNET (ML. LUU) ■ All • PI • ,.E*4 / 599. 

590 AIII (K) s All • '^SC (K) 

Go Tu (395. b35). IURDE 

59b ASS16- 340 TO N6HITE 
600 00 bOi; LOOK > 1 * 7b 

- IF (At I - BbODY (LOOK) )6l5»~ blO*~ 605 

60!^ CONTI mUE 

tEOdb ■ 100. - 

“So TO'630 '■ 

610 TEQ8B « 5. • (LOOK - I) - 27;. 

60 TO 630 
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61b if 62U» 625 - ‘ 

62a TEO«a ■ - 270. 

Go TO 6i0 

625 TEQHB ■ b, • ILOOK - 2) - 27 ,♦ 5, - SBOOy (LOOK 

ISbOUV (LOOK) - ObOUV (LOOK • D) 

630 TATM » To (M» UUU) - TtwHB 

PwInT 1516* iHtiA (K)t Mil* xCX. Xux* T^OBfrf-TATH 
1516 format (5HX*F9.2*t:i»*<»*F 12*4'E13.3»2F8 ,i ) 

60 TO NdRiTE* (660. 2fS) 

635 RRINT I5in» T*^IA (K>*'An» xCX» XUX" 

NUINE« ■ NLINF.S'* 1 

iSld FOhmAt (5HX*FV.2*Ei5.6tf 12,4^E13.3> 

***'6%0 K«K»1* ■ “ “ 

IF (K, I E. MANGLES) (»0 TO 2«' ' ^ ^ 

* — ■ ■ — ENO OF ANGULAR LOOP ^ 


GO TO (645. T25). 

— 645 ■ **K0 ■ mkn - - - 

SUM? s AIII (1) ♦ Alll (2) ♦ AIII (3) 

MKOLUr. » mKO ♦ LUOU 

print 1520. (MKOLUO) . TO -MKOLUO)lT)nMKOtTJO) • OOtttrtMKOtUO) 

ISLM7 . 

1520 format (Fft.l.Fli.l*FlJ.6.El2^*El^6/^ 

nlINEc « KLINES •“! 

IF (N| INES *LT. 3u) GO TO 65 • 

MlinE*; *0 ' __ 

UPAGF * npaGE • I 
LAlL ^lading (LUO) 

6S0 fiK«MK,LtVELS 

IF(MA.Lt.MKL) GO TO 281 
750 COnTI.iUc. 

NPAGE * 1 

NLINEc » 0 
65b COnTI.UE 

74b CONTI.Ut 

■* IF (UFLAG (18) 0) 60"T'''~^25 ' 

MKK » MKF *1 

IF (MkF .lT. 1) GO TO 770 

■ “ bO 66" Mk » 1 . RKF 




• 4 


MKP ■ MKK - mk 

660 VNET (MK) ■ FNEi (MK. 1) - FnET (MKP, 2) 

770 COnTI.UE ■ ‘ ■ ' •* 

PRINT lb?? 

1522 format (•iPRtSSUKE t‘NETJ_ A^M. TEMP. CHANGE*/) 

COnAT,. « 5880 . 

PRINT 1524, PH (1. 2). VNET (1), ZERO' 

1524 Format (Ffi. 1 .F U.4,F 1 o. l ) ..... . 

Punch 1526* PR (i. 2). TO (1. 2). FNET (W TT . FN£T (HK^? 2). VNET 
I (1), ZERO 

IF (MvF .lT. 2) GO T0_775 V 

bo 66*~ Mk » 2, MKF 

ATMCM « (VNET (MK - 1) - VNET (MK) ) / (PR (MK - 1, ?) - P« (MK, 2) 

1) * CkNaTm * 

print 1524. PR (MK, 2). VNET (MK), aTMCh 

MkP ■ MKK - mk 

Punch 1526, PR (MK, 2) . TO ( <K, 2^* FNET |MK, 1), FnET (MKP, 2 ), V 
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INfT «MR>t aTmCH 

lS2t> FOHMAt (lX*F7.1*ttX»Flt,lt3Fe.4tFa.l) 
66b CUNTIfUt: 

77i CONTIrUfc _ 

do TO 725 


• — ■-■- bi.ack body calibration loop ■ 

67o ASSIGN. o«0 TO Is!t 

- NL1NE<“« 0“^ * 

NPAOC ■ I 
temp i II 

AsS16». 710 to IOLACK 

67b 60 TO UT, (060* 665) 

660 AS5Ib*> ott'a TU 16T 

■ Call mEaoINO (3»* 

PRINT 1028 

1526 Format (• TtMP.**6X**iHHA01.»NC£«.8X.*RA0lANCE«»6X»«8A0lANCEV 

• I'lSAt^O/SO CM*.6Xt«*t/SW CM SR**5X*»(N0RMAL)*/) ' - - 

665 SUM6 r Sums ■ 0* 

IF (Nm2 ,lT. NUI) 60 TO 780 •• 

00 70« r« NUlt‘NU2»* NOJ : 

UNU « I 

P6R3 ; 0 • UNj / (TEMP ♦ 273.16) 

A1 M i, - --- 

IF (P^Hj ,lT. 6V0.) AI ■ p • UNU • • 3 / (PI • (EXP (P6R3) - 1.)) 

■*iw « T • .1 ’ 

- Iw « wW ■ " 

FLTRItI ■ FLIk dw ♦ 1) ♦ (Ml* - IW) • (plTR (IM * 2) - FLTR (Iw ♦ 

11 ) ) 

‘ SUM66 » Al * FLIKIII • UNU 

SUM55 » FLTHIir 

69o IF ((T ‘EO* )^U1) •UP* (I .E(2. NU2))695, 700 
69a * SUM66 ■ SUM66 • »5 “ ' 

SUMSS « SUM55 • lb 

70U SUM6 * 5Uh 6 * SUi-io6 

SUMS ; SUMS ♦ SUMSS ■ 

7ob COnTIkiUc 

780 continue •• 

GO TO IULACK* (/lO* 235* 255) 

710 SUMN 4 bUM6 / SUMb 

SUm7 i bUM6 * CC22 

print 1530* TemPi SUM7* SUH6,~SWN 

1530 format (F7.liEl«.6»2Fi4,7) 

NlINcc » NLINES ♦ 1 

- IF (Ni INES .6E. b0)715» 720 “ 

71b* NL1NE« ■ 0 ‘ 

NPAGE • NPAGE * 1 

ASSIG»i o80 TO 1ST - “• ~ 

720 temp , TEmP ♦ TU* 

If (TfMP ,LE* TT)675* 725 

725 WetURm ' " ■* 

End 
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Sur-HOnTlNF TahPHnT <AL* IhOAi ) 

CUmho» /TRANS> /Ii«IA1» ImaO* Ny^At TT (50)» TCOlt TCOO* NTCO» TTC 17 

l4) tTOTl’TOZOtNTOZ.TTOZiHl) 

UlMENclON AL IMUAL“(2» 

UlMENclUN IFhu Is) ♦ VAI, (s) 

ENTRY VLNRRNi 

PRINT ISOO* IHOAL ■" ' ■ ' 

1500 format UmI * J 6Xt2A«//2X*5(» FREQ LN(l>*)^) 

UO 10« II ■ I* 


100 

1502 

105 


110 

115 


1504 


i 

I 

‘ 1506 

120 

I 

I 

J • “ ■■ 


150u 

12s 


" 13o 


1510 


KK . tI “ 

UO lOn JJ ■ It 5 

iFHQ iJJ) « (XK - 1) * 10 

VAE ( iJ) B AL (KK) 

KK B vK * 50 

PRINT lS02t (IFHQ tilt VAL (T)t I b It K ) 

format (2Xt5 ( 1 7 tF7,2) ) 

COnTImUE 

00 lU II » JOt i>0 

KK B T I 

00 11- JJ « At ♦ 

IFRQ (JJ) » (KK - 1 ) * 1 0 

VmL ( iJ) * AL (KK) 

KK B kK t 50 

PRINT lso2t (IFRQ (D* VA L (T) » 1 ■ 1» A> 

CunTImOE 
Rf TURf 

Entry t«ater 

FURMAT^^lMl*3lXt*RATER VAP0Rtt//7Xt*LN(L) TAU*t2 (9Xt*LN(L). 

1TAu*)/> ^ . 

Do 12o I B IBt dO 
13 B T ♦ 20 

vlna , I * .1 - J.a — 

VLNd ; vLNA‘* 2i 

PRINT lSC6t VLNAt TT (l)t VLNbt TT (13) 

format (21Xt2(Fll.lj^l0.4n 

ConTU'UE 

00 12e 1 B It 10 

12 B T ♦ 20 ‘ : 

13 B f2 . 20 “ 

Vlna i i • .1 - 3. a 

VLNH - ylna"* 2. " 

VLNC I VlnB ♦ 2. ■ ■ 

PRINT IsORt VLNAt TT (Dt VLNBt TT (12)# VLNCt TT (I3) 

format (3(Fll.ltHoi4)) 

COnTU.OE 

00 13 : 1 e lit 17 

12 ■ T ♦ 20 — 

VLNA i I • .1 - 3.d 
VLNB , vlna'* 2» 

PRINT IsOHt VLNAt TT (I) t VLNBt TT ^l2) ■ 

CONTlfUt 

HtTuR* 

Entry tcarbon 

PRINT Iblo 

format (lHlt2QXt*CARB0N DlOX»OE»// 

1~ 7X t*LN (L) •*5Xt *TAU*i 2 (bXt *i^(L) 
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SU0HOiiTINE HtAO"AT <LUO> ' ~ " ' 

print HtAOiNG A'lO INPUT parameters 

COmmOi" /head /lilAi MOaIE* N ATEt WLl» T^LZ# IF1» IF2» FNO^i Fno 2» 
1 NPACrt OMJ* NPfWf IOHIIE --- 

C 

DImENcIOn NMtAD (3> 

^ - 
r 

Data ^nheaD « h » upaaku, sh >ownwaro* bhUB ^^y^***_ _ _ 

PRINT l3(»0« ISTArMOATt. NOA^E, NParE 

1500 format ( Ihl ♦•STAJlUNi Aflf 5 < t*OATE * •» A0*b6X««CAi.CUtATEO# •»Afl» 

l'5Xt«oAoE ♦»13) ” ^ 

~ print lS02* oiLl* 1 aL 2« IFli TF2 " 

1502 Format (•OMAVE UENoTM region from •,FB,t>»* to •»A8» 

l"* MlrRUNS WITh'fIlTLR ■ *^2 b)_ ^ 

LL » I U" 

PRINT 1304, FMU1» FNU2, ONU, NMEAO (LL), NPTR 
150a format (TXt^^AVt NUMBERS FROM *«Fd,2,*- TO 4» FB»2,* BT *,Fd »2»2 8Xt 
■■ i’mh,* RmOIANT Pu-Ek FqW LAYERS*) 

PRINT ISOb 

150o Format (1H0*11X»** 0**lSX,*RAOlANCE*»llX»*RArER VaP0R*»9X» 

■ r •OlFFtRENCE*/12X',3MPPM,i4x.»W/SO“CM~SR*»30Xt*W/SO CM SR*/) 

return 

END 






OO 


subhoiiTine MtAOiNG (uuy) 


■ PRINT heading and INPUT PARAMETERS 


Common /head /I^TA, MDAIE, N^ATEt WL1» IF1» IFZt FNUU fni^» 

■rNPAGrt ONU» NPfiiribHiiE 

piMENelUN NHtAD (3) 


DATA (NMEAO • BH uP*‘AHi>» 8 H'.OwNhIARd, RhBB CAI.IB) 

print IbOOt ISTAtMOAItt NDArti NPAGE 

ISOO Format (1H1 »*STATIUNt AtAgfbvt'AUATEt •# A8t66X»«CALCULATE0» •»A8» 

1 bX*»oAGE *tl3) __ _ __ 

PRINT IbOPt'ALl* IWL^* IFlt »F2 

1502 format «»0WAVE length REGION FROM *,F8,9»* TO *tA8» 

!"• MIrHONS WITH FltltR*" •»2 -b) 

LL * I UD “ 

PRINT l30Af FNUlt FNU2t UNU» NHEAD (LL) . NPTR 
1504 format (7X*A"aVL NUMbEHS FROM ••F0.2«* TO *»FB.2t« BT *tF8.2 *28X * 
— 1'A8»* RAOIANT pOikEr COH *,13.« LAYERS*) 

GO TO (iOO* io5)* IdRllL 


lOU 

1506 
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1508 


Format (*0PRtSSORt*i9Af *TFMP;*»5 X*amIX- pATl0*,5Xi*0ELU*t5X**lRRA0l 
1 ANCE*. SX ** angle* •3Xt*RAUlA'.'<CE*»8Xt*C0 2i«8Xt«H20*t7Xt •BRIGHT ATm 
2*/4yx.*«/S0 LM** l5X»**</90 j«M $R** 1 A»2 <AX****M^SO Cm*) ♦ 5 A»*TEMP*tAX, 

3 •C0Ro») ~~ "■ ' ■ 

RETURn 
print 190P 

format («»0PRtSSURE*f9X*«TEMP"*i5Xt*MlX'oATl0«,5Xt»DELU*.5X»«lRRA0l 

lANCE*.5Xt*ANGL£**'»^»*”^^‘^l^'^C-**®^**C0 2»*»X**H20* 

2 /a9X,*m/SQ Cm**' l5Xt*i*/9Q CM sB*» lAf AXJ*£^''SQ Cm*) ) 

HETURm" / ' '■ - - - _ — _ 


onnonooo 


FUNCTtON nSUoS I*** VH, Np) 


SEA«CHt5» HHESSURt AHHAY Tq DETERMINE SUBSCRIPT FOR “ 
spttiutoVRFSs Re 
p tS the Specified ph£sSu-e 

VP IS THE'aRHAY of PheSSUpE VALUES 

NP IS THE NUHijEK 0F*HR£SS .he VALUES 
NSiiBS IS The UtSlKEU SUBSsHIPT 



UlMENcIUN VP (IV) 

If ,^p 6 o TQ-utr •* 

Uo lOr 1 s i* NH 

IF (AfiS (p - Vp'(D) .LT. ,ui) 60 TO lo^ 

TOO continue ~ ' ~ ■■ 

llO CONTl».Ufc 
NSUHS « 1 

RETURm 

los NSUBS > I 
HETUHm 
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FUNCTtON TAUw (VnUJ 


*ATtH Th^nSmiTTanCc FOH NU m (2 5*797.5) BY 5* 

COMHOfc^ /XYZ*/ /X* Yt Z* C (9t X20) 


■w (I) 

in ♦ 

2 * 

WS 
IL 
*tL 
WA 


C 


(C 

<U( 


(4,* 

I) 


I) ♦ 

• X) 


C (I* 1> ♦ 

♦ y * (C tfc* n 

C |b. !)■ • X) 

■ |VNO ^ 202^5 )■ • .2 
« wS 

■ M (iL * n 

s -;;L' '*~n<5-jru.r T au~*“2) 

TAUW ■ txp ("- txp (Wa)) 

AETURm 

end 


r • 

♦ C (7* I) 


• (C U. I) 


“WD” 


[| 

il 

0 



FUNCTtO.N TAUCO tVNU) 

CAHUUN OtOAIuE iMANbHiTTANCE FO NU • <507.S*8i7.S) 8Y S. 

COmhOm /XYZ /At Y, zi CCO <9, 63) 

W (I) » CCO «lt I) ♦ Z‘» (CC ) |4t I)-* X • -(CCO t5» X)~* CCO (9t~ 1“ 
1) • Z») ♦ ecu (Jt I) • Y ♦ X • (CCO (2. 1) ♦ CCO (5. I) • Y • X • 
2(CCO |7t I) '* ecu (bt I) • Y) ) 

- WS » (VNU - SOTIs) • .2 - ~ 

IL « wS 

bL ■ t. (IL * 1) 

WA • U.L'* (WS---IL)“^ (2 (IL -♦“2r -lICl 
TAuCO ■ tXP i - tXH ^«|A)) 

AETURm * 

tNO " " 


ORIGINAL PAGE IS 
OF POOR QUAUTY 


'jUawOi.I INE OUfiMU 

COmmc). /XYZ*' Ci^ *9t 12.Y) 

Data ,c* = b.i2i«* .9i<*<:t .v ii. - ,437 a* 

1.), 6 W'i, 1.02»'* 0.'. .ObAoo, .n<#a9 

^ . .t. 1 ^ . li \ . ti - n»^ii m . - . nR4 1 . . n976b« 


,049b6f 0.* .Ol«lA* ? <0 
9i ,02A4* i (0.)* 5.6094 
I. 2 (0.) . A,tt9Si>. .9SA. 
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RRPRODUCIBIUTY OF THE 
URUilNAL PAGE IS POC« 




SUBHOiiTINF 
COMMOr 





ATS, ,75<iW» 0.* O.uli'* 

bi.l2bA« - ,0/o7» T.d’^aOt .0* 2l» 0.) 

DATA, (CCO(I) 

w ■ J.» ,S0f.40N - - 6.929«» «•’ 

112, 0 . .A47o, .02279, ".01200, - .1162, - ^^^296, .9*j7?, 0.» 
221, ,r3299, 0.,"0., .OJOISB, - .06665, - 5,65, 1.0686, •028l* 

319, .«3'*a5, 0., 



A3<»8, 0306, - ,ww— .... «.. - .u"---’ 

- ,oo30i»-i“,000057) 
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-.4 


bUHHOi.TlNE OUMMT* 

COmmUm /XYZ /uX* UY» \il* ecu <9t 63) 
Uatai iCCUlI) tI>^V7iJ6/) 


U.t " IMUIDJ. - •UUlUC'<t «•’ - 

6» l,5^6*t ,0J761t - .G0Y6t - .00634, - .002023. - .00233. - 1.41U6 
7. .75?. .3614. 2.2 o 72. .049*.'. - .0049, - .01173. - .0031,04. .003o 
bl. - T.96?3. .8619. .J2/. 2.H036. .06731 • ” .2118. - .01*192. - «®0 _ 
94801, .01569. - 2.546. 1879. .2863. 2.9)8. .07776. - .2568, - .021 
*24. - .005440, .02821. - 1*6 3. .3043. ,2803» 2.9613* .06098. - .2 

0836. _ .03415, - .005813. .0)779) _ _ 

U*TA((CCU<I»*1«J88,479) 


8. 0.. ” ,04328. - 0.7219, .9^8. .0404. 

9. .00/’-046. - ,0«89. - 0.0114. .9978, 0, 
*05637. - ,075o 5^ - 7.4221. 1,0581. o.. 
84869.*- .00853. - 7.8763, ,8^91) 

Liar. I .n iMT) .Ts<»HU<OOn 


B.C470. .04*1. u.« V.. .n 

.1331. .04566. 0.* 0.* .00 

5, .6383. .04162. .012539, 
879. 0.. 0., .O 38 I 7 , ,0140 


6,0192.. U.. " .'?i06* ” o^o-.3. i.uuoo. - 
”718311, .00973. - .000534. o.T " 7,7652, _ ,5169 , 3 (0, 
8b. ,0T1 o44. • .28153* - 8 . 31 ^ 8 . .6049* o»* " .8062* 
995 , .»5332< .018 866. - .3415, - 9. 4568 , 1.Q304. 0. * 
A. .0l(.lo4. 0. ) "’ 

END 


— SUHMOtiTiNE OUMMY5 — 

COmhO-. /ALFLTH /*L (229), FLtH (229) 

Data »AU « llio* T.. 40 . l.*2. 1.90. 2.04, 2.20» 2.30. 2.40. 2.47. 2 

ll55. 5 . 00 . 2.65* 2.67. 2 . 68 , 2.68, 2.69, 2.69, 2.67, 2,65. 2.63. 2 

2 * 61 . ?.o7. 2.53* 2 , 48 . 2,4b. 2*40. 2.35. 2.30. 2.23. 2.15. 2.10. 2 

37o 3, T.96,- 1.07, 1.80, 1.72, i^bb, 1.57, l.SO, 1.42, 1.30* 1.20* 1_ 

~ 4 I 12 , T.03. . 95 . .U7.'.77, . 67 * .60, ,52. *^2. ,33. .23. .15, .05. 

5- ,03, - .15. - ,21. - .30* 1 ,42* » .55. - .65. - .74. - .85. - , 
697, - 1.07, - 1.20. - 1.31* - 1.42, - 1.55, - 1.66, - 1.79, 37 - 

■7*8.00)) 

UAy*( (*L(I) *1 b1IU.218) 

x" * 17 ( - 8.00)* - .65* - .50* - .26* - .13* ,05. .2 

~ 10, ,3c* , 50 . ,63. . 73 . . 05 . ,97. 1.07. 1.17. 1.27, 1.37, 1.45, 1.5" 

25, 1,*.4, 1.75, 1,84. 1.95, 2,05. 2,15. ?.28. 2.40. 2.46, 2.42. 2.2 

38, 2.T6, 2.01, 1.75, 1157, 1.48, I. 56 , ).87, 1.83, 1.98, 2.08, 2.1 

— 4b, 2.52. 2.2». 2.29, 2r28, 2*26, 2.23, 5.17* 2.12, 2.08, 1,97. 1.0* 

58 . 1 , 70 . 1 , 80 . i.5o* 1.38. 1.27. 1.14* ).03. ,91. ,02. .71. .61, , 

650, .*1. ,31. ,c2. .14; ,o7. - .02, - .)2* - ,21. - .28, - .40, - 

7 . 49 , _ .58, - ,66. - ,7b. " .B6.‘- ,96,“'- 1.07. - 1.17. - 1.28,-“ 
81.40, - 1.50, - 1.03, - 1.76. - 1.80, - 2 . 06 , - 2 . 14 , - 2.29» 

UATA( (ALd) *1 >217,229) 

X « 2,42»- --2*Tj5." -"2. '^2.'”" 2.8T. -■3,03, “ 3;20* - 

1 3.38. - 3,55. - 3.6?; - 3 . 85 , - 4.03* - ‘♦.22, - 4 . 4 I) 

UATA ,FlT« > 22? (l.O) ) 

— END — ' 



ORIGINAL PAGE IS 
OF POOR QUALmr 
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SUHHOiiTiNf; OJHMT^ 

COmmOm /Ai cal-/ ALC(^^V» .AL•l',^^9^ .AlOZ«‘* 29) tUZHXMn9) tOZP»<Sll9) 



C iOU(L> fABLbi* KUH CO?t mlNUUN, AND O/Oi^E (FMUH EuSAS&tP) 


OATA 

,AlC ■ 

36 (0 

• 1 • 

• 

97, 

- 

4..I5 

1 , - 3 .J 0 , 

• 

2 . 

63. 

• 

2 . 06 , » 

1 , 

5 

10 , - 

T>'j5, 

« ,61 . 

• • 

21 . 

• is: 

. • 

18, 

.43, ,32. 

0 . 

f 

• • 

34. 

- .72, 

• 

1 

2.08, 

. 1.50 

* • 1 • 

68 . 

- 2 . 

2 V, 


2.67 

, - 3,08. 

m 

3. 

52. 

• 

3.93, - 

4. 

3 

39, - 

4.63, 

14? (U 

.) ) 













OATA 

*,ALw ■ 

68 (0 

• ) t 

• .81. . 

• 

71, 

- ,73. - . 

.75. 

" . 

76. 

- . 77 , 

w * 


178, 

- .60* 

- .81. 

• • 

bJ. 

• • 

>4 , 

• • 

86 . . .87i 

t • 

• 

89. 

- 

.91, - . 

92 

f 

2 - . 

9*. - . 

93, - 

.97. 

* • 

99, 

- 

1.01 

* - 1.02. 

- 

1. 

04. 

. 

1.06, - 

1* 

0 

37, - 

” » U8 , 

” 1.10 

• • 

1,11 

f • 

1. 

12, 

- I.l4, - 

1. 

15 

» • 

1. 

16. - 1. 

17 

f ** 

4 - 1 

. i 7, - 

1.18. 

- 1. 

18. 

- 1, 

.20 

f • 

1.2«, 4 ( 

• 

1. 

21) 

t • 

1.22* 5 

1 ( 


5- 1. 

2.), - 

1 ,20. 

- 1. 

19, 

- 1 

18 

9 • 

1.16, - 1. 

.14 

• 

108 

(0 

.)) 



OATA 

(AL07« 

97 (u. 

) ,-2 

** 

I .4, 

1 * • 

8*>.4 *’•,‘>5, ,2i 

>.3 

9 • 

5.- 

*1. 

,5. , 35 ,“ 

“■ ■■ 


♦-.7, 

" V » 23 , - 

1.3,-i 

#««# 

1 • 

•••• 

1,3. 

•••; 

,-2 

>• 

.5* 

114( '.)) 

f 

• • 


• •• 


• • 

• 


c i.iaing hatio vs.‘ phessu f Tables fo« ozone ■ — ^ 



Oata (OZPRS>iuS 9 . »lUi\l.'* 6 l 5 .*'e 06 .« 356 .t I 9 <»* tlA 5 .» 103 < tS 5 .i 29 .*t 

OATA rOZMXN>U. 0 *i* 9 * 7 .<»tH.A; 9 . 0 f lS.,lB.,A 3 .* 300 «t 733 ** 1037 .t 
• 1154 . # 1180 . ♦llWJ..« 97 .t 702 .. 500 . t 370 .t« 9 .> 

End - --- “ - ^ - - - - - 


- SUBWOiiTINE C)yMM)r 7 - 

COmmOm /TftANS /i-AIt T-mD. Nf-A. TT (Sfl)* TCOl, TCOO# NJCO. TTC <7 
14) ,T07l.TOZOtNTOZ#TTu2(*#n 

* 

c tAU TAdLtS FOh vaPOh# C 02 # ANn OZO'^E (FROM ElSaSSER) 

* 

OATA (T-AJ * - J,7)» (I-AO ■ .1). (nT*(A » 50) ■“ ■'* 

OAtA ,TI « 1., .99V. .9967. .9932. ,98R7. .9833, .977?. .9705. .96 

133 . . 0536 . . 9 ^ 73 . . 9 J 89 , . 92 no# , 909 a. . 8978 , . 8 R 51 , . 8711 , 

2 ‘,« 55 a, .h 3 R 4 , . 8197 , ./ 9 R 4 , . 7753 , . 75 ?» 0 . . 7225 , . 6928 , . 6610 , . 6 " 
3272 , . 591 b, . 6341 , . 3132 , . 4 -» 5 o, . 4 -^l 8 , . 3919 , . 3497 , , 30 T 6 , .2661 

4, . 22 <: 8 , .1877, .1316, .1198, .(j920, .0s*<7, .0500, .0357, ' .0250, . 

50168, .ylo3, .0079, .0) ■ - - 

OATA jTCOI ■ - 3. 2), "(Icon ■ .1), (lyTCO » 74) 

OATA ,TTC ■ 1., .9997, ‘.9991, .9982, ,9o70, .9955, .9937, .9916, . 
19892, .Vb 6 S. .9835, ,9802, . ,766. ,9T?7. .9685. .9639, .959, ,9538~ 
2. .9402. .9422. .9333, .9?9, .9218, ,-»l 4 i, ,9o59, .8972, .8879, .8 
378, . 06 / 7 , .6561, . 677 , .831, .817, ,Ho)9, ,7850, .768, .7491, ,72 

488, .70/1, .6839, .6392, . 633 , .6053, .9762. .5458. .5142. ,4815. 

5 . 4478 . . 4133 . ,J/ 6 <:, .J 7 ;> 8 , .3075, , 2727 . .2389, ,2o*>6. .1763, .U 

667 , ,T2J, .lUu2' . 0801 , . 092 ^, jOAS, , 076 . ,0266, . 0196 , * 0 lA 7 , .0 

7113, .0069, ,uO/l. .0030, •0*42, ,0028,~«00l4, ,0) ' 

OATA (TuZl»-4,3l . (TOZO«.l) . ( T0Z>41) 

OATA ,T I07« 1 . . ;99<i, ,967, ,97t, .968, .96, .952# .944, .936. *92 7, .917, 

'.906,.897,.88l,.Eb6..8H9,.83,.809,--7d6,.761,.734..705,.673, * 

•.683, .6, . 559 , .313, . 46 b, .418,.-)66,.312,.257,.202..1S1.,108, 

•,U74, 048, .029* ,0l3, .003,0. 0i 

End ' *■“ 


EEPRODUCIBILnY OF THE 
ORIGINAL PAGE IS POOP 
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SUbWOiiHnE *»*<OFlUfc{rt»NHlK,AL '»NU3) 

COmmO»> /PA«AMS no («»)• AWlN. XMAX, YMIh. YMAX. NC^C 

COMMO»/HT»./AH(2yo».xrifJO»»x (2oo» .jACznO) 

C0MM0. /HEANS/PBAH(lbNiaAR(l5» tlHEArMB) 

COmmO.. /PKOF /rtAOC (lal* MAO'* (I5»t NFBOL <15)» NF«QJ U5)# NFHQ 
blMENclON’XX (l&)» UX (iS)« F.E (15) » Xpr (l^)» TAE (15)» AXE ( 15> 

dimension KX (lU), PBAK (15). T8AH-(l5) 

DlMENclUN XTh (ibt 15) 

A»357.9U0835 

16U KX (l» • I ' 

KX (2) ■ 1 

KX (3) ■ 1 

“ KX (A) ■ 

KX (5» ■ 10 

KX (6) ■ 500 - 

REau T53A, li) 

15JA FOPHAr(aAlO) 

HEAD ibJA, I«EAU _ 

PRINT l3Jh» ImE**!]^ 

iS36 format 1 1h 1 taoXtBAlo/) 

MEAD T5JR. XMlNt XWA X. Y MIN. TMA X 

l53b format (irFB.2) ' 

REAU 75-*0. NCYC 

15Av) format (ID _ _ 

call (A. a* Ct y) 

ME > ^ 

HEAO T5«*0. MtX _ 

iF (MeX ,EO.*0> 00 TO JSS" 

ME « T 


~ iSAi Format i/t pRtS. TE-^T*, w EXAcTk/) •• 

36 S read TS**A. XPE (Mt). TAfc (ME»» AWE <mE) 

I 5 a*» format ( 3 F 8 . 2 ) _ . 

If (XrE (MED 373 . 375 . 37 o 
370 tAE (mE) « TAE (Mt) * 273.16 

Print 13A8» ape (me)» (ae (Mg), awe (HE) 

ME • >'E « 1 
GO TO 3 o 5 

375 me » . E - I 

• MEl » Mt - I ■ ■ 

IF (MfI ,lT. D 00 TO A 75 

bo 3 M« 1 * It Mtl .. 

- PHaR (I> ■ .5 •'(XPE’d) ♦'APE (I * D) 

380 trtAR rD ■ * (TaE (D ♦ TaE (I ♦ D) 

475 continue ■ 

Call tRCOT (APEt TAEt MEt D 

1546 Format </X PRtS. ___ 1 _ initi al GuESSp/) •• 

IF (NnT- ,LT. 1 > GO 10 480 

bO 3 V»' I « I* NPTi* 

PRINT 1548 t AP (i)t TA (I)i_xW^ (I) • 

1548 format (/F 9 .l.Fi 2 . 2 »PHiA) 

390 COnTImUE ■ ' 

480 CONTI.Ut 

print l 5 S 0 t ALO 

1550 format (//• TnEIA ■ ••F5,1) 

C (PROCESS) 
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1556 

349 

1556 

4UU 

465 


405 

490 


410 


415 


HHint 165? 

POMMAt <//?« I»jX»mtUUtNCie^*»2X*R4kOU»'CE*/) 

MtAO S59A, NPHQt HUi 
POHHA-r n?tl4) 

IF (NrW'^ .Ll* !• ^0 •♦•*5 

l)U AOf I ■ It UPHU 

HEAD ?556» NrHOL (I)t NPHOU (lit RAOO (!> 

FOHMAt (2l5tFl5.4) 

PHlM I550t it NFHUL *l)t NF (UU (I)« RahO (I> 

FrtftMAr (/12t2I8tri4.91--- 

CONTIfUt 

CONTIt-Ot 

-N ■ M -“1 

IF «N .LT. 2) GO TO ^90 

UO 40r 1 ■ 2t N 

XX (I -'!) • TA tlr 

ux (1 - 1) ■ lOt 

CONTUiUt 

- V , — 

KX (7) ■ KX (8» ■ XX (9) 

NITER ■ 0 

• z - 1.- ■; 

CAUL rALC (XX, ttt V< 

ASSlGf 410 TO KHRT 
Gt) TO 4 15 ' 

ASSIGm 450 TO KHHT 
00 45k nITEH » It 2 u 

CAUL kINMYZ (XXt OXt N -"It'^Et NFRQ , 2 » • “Vt~ XTH) 

PHINT l660t V, (XX (J) t J 4 7t 10), nItf*^ 


Kv (10) 


P*,4,X,4 T U ■ •,E9.3t4l4,10X,*NlTER ■ •tl?) 


1 

IF (N ,LT. 2) GO TP 495 
00 42-t I » 2t N 


420 

495 

TA (1, ■ XX (T - 1) ”■ 
COnTImUc. 

ir IM .LT. 1) GO TO 600 

• • 

• • 

425 

t)0 43 k I « It M 

IF ((T .EO, H) .OH, (I .EO, 1))425, 43o 
Print i 6 t> 2 , xp Tn , ta (d 


1562' 

format (/F8.2,2P9.2tFl4.9) 


43o 

<»0 TO 4.i5 

PHTNT l96?t XP (1), TA (I), oX (I), EE (I) 


- 433- 
500 

CONTI. UE ' ■ 

CONTI.'UE 

PRINT 1661 _ 

• • 

■ • 1561 

format (//• PSAK -- it»AR*J 

Mm s V *• 1 

IF (M.. .LT. 1) 00 TO 605 

«• 


bo 44 '*■ 1 E 1 * MH 

T6aR (1) * .6 * (Ta (1) ♦ TA (I ♦ 1)) 
PRINT I662t PHAK (I)t fbAR (,) 


"440 

505 

CONTINUE 

CONTIfUc. 

Call tPlot (xp, ta, nru, z ) 

• • 

" 445” 
450 

GO TO XPHTt (4l0, 450) 

IF ((.HER .GT.*2) .Ano. (V .LE. l.E-9)) 6° TO 460 
IF (Xv (10)*»E0«. U) 00 To 46 


-455 ” 
460 

■C0NTl».Ut 

call pHIPlT (Ot Ot Ot Ot 20) 

RETUHm ' 
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subroutine CalC (A, t. U) 

COMMUt / hTk/AR<2u*J) tAK^JO) «A J(200) tTATP'^O) 

COmmOm /HkAO /ISfA, MOAIE, NfNAJEt WLl* J»<l-2» IFl» IF?« FNUl > FN L>2> 

I'NPAGr* ONU* NPIi«t IdHlIE " o « 

COmmOu /PAHAflT ^M» NANucESf wUl» NU2» Nu3» MKt LUOO* NCOS»l» r» Q* 

lPI«LUn.N02SW ' _ - _ 

COmMOm /WATEHp /KNObS. FZW. ^2T, VLAmB " 

COmmOm /prof /HAOC (IbTi RaOO IlSTt nFROL U5)i NFRQU (15>* ^FRQ 
COmmOm /VALUt /All 


UlMENclUN X (15) t t UP) 


N ■ M - 1 

IF (N .LT, 2) O TO ilW 

bo 10« I » 2t ) 

T» (I, « X U -'IT 

AT (I, ■ TA (I) - 27J.IO 

100 CONTlfUt _ 

“ 110' C0Nfl».u£ 

0 » 0 . 

ONU * NU3 

IF (nfR'J .lt. i)*GO To-iiy 

00 10c NF « ii NFHU 

NUl » NFHqL (nF) 

NO? * NFHOU (nF) 

KnomS » RAOO (NF) 

Call rAtwAT (xif_Fi«2t i) 
RaOVAI ■ FWZ - - 

RAOC (NF) » All 

E (NF> ■ RADVAL 

U 3 U ♦ RAOVAl 2 
105 COnTIhUL 
115 CONTI^.Ut 

"Us SrRT NFHQ) 

HETUHm 

£no _ 


J 


age is 

quality 
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I 


lOU 
— “*S« 


I OS 
ilu 


lis 


120 
460 

12S 
130 


13S 

Uu 

US 

463 


ISO 

15b 


J 


-^SUhMOmTINE M1NMYZ-(Y» UX» NN, E» MMt 4, KKt Uf-XlM)— - 
dimension X «15>» 0* <ib). E »15)* EE 05) » XX (lO)* OE (l5)* DDE 
l(l5)t XfM (16« 13)1 ATM 05* 15) » XTNA (15* 16)* OEL (15)* Y 05^ 
Closed (a* x* b) ■ amIni (amaxi (amini (A* b)* x)« amaxi (a* on 

VN ■ r 
M > Mm 
N a 
V a U 
Z a W 

IF (N ,LT. 1) GO TO 450 
00 lOf* L a 1* N 

X (L) ■ V (L) _ 

COnTUUE ~ 

2o ■ 7 

VO a V 

■ NUM a 0 ' ' 

KX6 s Ka (1) 

XX 17 ) a 0 

KX (Bi* ■*() 

KX (V» » 1 

KA (1,) X 0 

“PEH a 0.1 • Ft 0 AT“(KX Qll 
ZM a ^.JOl • FLOAT (KX (5)) 

2 a OaER <Zm> U)_ 

FM a n 
ULM a KX (4) 

IF (O|M)110. 105* 110 

DlM a IjO. 

IF (N .LT. 1) GU To 433 

DO 19- I a 1* N 

IF (l)y (I))115» l9o* 115 
NOM a NUM * 1 

X (I) a X (1) ♦ OX (1) 

Call rALC (X* Et* VN) 

0 a 0 

IF (m'.LT. 1) GO TO 4oO _ _ 

"do 12- L a 1* M 

ATM (| • NliM) a (EE (L) - E Cl.)> / OX <!) 

0 a 0 ♦ XTM (L* NUM) * • 2 

CONTImUE 

If (0)U5* Ub. 135 

X (1) a X (I) - OX (I) 

□X (I, a 0. 

NUM a NUM - 1 

GO TO IVO 

IF (V* - V)l40* 153*" IsS 

IF (M ,LT. 1) GO To 463 

do 14^ L a 1* M 

E (L) » EE <E) ‘ 

COnTUUE 

KX (B> » 1 

bet a 130. * (V - VN)-/V ' 

V a V* 

IF (drT - PEB)1V0* J30* 150 

Kx"(7) » 1“ 

GO TO UO 

OX (I» a - 


• • 


• • 


• • 


• • 


• • 


• * 


• • 
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IF (Ul M» 160* !<»«• 


X (1) » A iV ♦ uX 
GO TO l'*n 

X (I) * X (1) ♦ GX (1) • ?. • ' 

CALL rXLC <X» Et» VN) 

D s 0 * 

"iF (M .LT. 1) GO TO ^JO - . 

bo 17- L « I* M 

XTm (| * NUM) « ifet <L) “ E (l.)) ^ OX 
D ■ D ♦ XTM (Lt’NUM) 

CONTiMUt 

If i'5» i«*o 

X (I) » X (I) ~ox m ~ 

GO TO IJfl 

I IF <V»- - V) l'»0« iGb* Igi* 

r X (i)-»-x (i)*- ox- tn 
Ox n, ■ ox ( 1 ) • .s 

I COnTImUE 

i -continue 

If <N|.M)195* I93t 2lU 
> Kx (9» ■ - 1 

- GO TO A20 “ 

0R^G1NAL MINMY uECK INSERtED MERE. 

; IF (M ,LT. 1» GO TO 

bo ? 0 c LL ■ I» 

t EE (Li ) ■ XTM (LL* NUM * 1) * E (LL) 

} Contu.ue ' 

(30 TC 2 70 

J If (M .LT. 1) GO TO 

00 22« LL “ ” 

3 EE (Li > - ABS (t (LU > 

j XTM (i L» NUM ♦ 1) « 1. 

j COnIImUE 

CALL cOMT (EE. M) ‘ 

Emx « Et (M) 

£MN « Et (M ♦ I - NUM) 

" IF (EmN .eQ. (i> go to 225 

IF (IrOJ. •"emN - EmX) 2J5» 22^* 225 
h IF (M .LT. 1) 60 TO 
~ bo 23.- LL “ I • ‘ ■ 

0 EE (Li ) » E «LL) 

3 CONTI.tUt I » 

- GO TO 2s5 

b EX s *.9077533 / AlOG (tMX / EMN) - 
IF (M .LT. 1) 60 ro ‘»90 


, 1) 00 TO '♦BO 

a I, M 
ABS (t (LU) 

siUM * I) m i. 


► 1 - NUM) 

(i> GO TO 225 


IF (M .LT. 
bo 25- LL 
If (ArtS (F. 
FAC *: *AOS 


1) 60 ri 

t l, M ■ 

(LL)) - 
(t (LL) 


£MN)2S0» 240 

/ EMN)) • • lix 


IF (N|.M .lT. 1) 60 TO 495 
bO 24t NK » 1. NUM 
XTM (| L» KlO ■ * 

COnT I k'UE 

XTM (| L« NUM ♦ 1) ■ FAC 
EE (Li ) » XTM aL. NUM ♦ 1 
"COnTUuE ' 

If (N..M .lT. l) GO TO 3J0 
bo 26e II « 1. NUM 
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## 




i 



J 



U (Ni.M ,lT. in GO 10 SOS 

do 2 bc J B lit nom 

xtn (t1* j) ■ 0; 

IF IM ,LT. n GO TO Sio 

do 26« KK B I, M 

2b0 xtn <tI* j) ■ X|N (II. J) (KK, Il» * 

310 CONlI.iUt ’ ■ 

2»jd xtn (.1. ID ■ XJ^N (II. J) 

50S CONUtxUt 

500” CONTI‘-UC ■ 

27y IF (NiiM ,lT. 1) GO TO 515 

do 27c LL “'I* *'*0M 

XTNA *(LL. NUM- ♦ -17 ■ TT. ~~ 

IF (M .LT. 1) GO TO S20 

275""XtnA^VLlI*NUM*^-1)~t«“X TNA (Lr*r-NOM--i Tr"<^”XTM -UT »“l. L ) 17 

520 CONTIrUt 

51b CONTl.<Ut — — 

If INhM .LT7~n”G0"TU-525“- 

do 2«- II b'i, ^UM 

• IF (Ni.H ,lT._1) GO T^S^O 

"do 2H' J B 1 . WOM 
2rio XTNA fll. J) * XTN (11. J) 

530 CONTIfUt 

52b CONTlMUt ■ " ^ 

Call r.AOSSEL (XINA, lb. NUH, NUM ♦ 1, LOL> 

_ If (LrL - NUM)2 o^._2'i>0. 285 

28b KX (9> B 0 

P8INT IbOO. NUM. LOL 

1500 format (A SIN6U UAW MA THIX. 0b0E RA.I2.F» 8 ANKA I2»»A) 

" Go TO a20 
290 NS » »“ 

295 IF (N .LT. 1) SO TO b3b 

do 30c K B 1. N 


• • 

• • 
• • 
«« 


#• 


• « 
#« 

• • 

• • 


#• 


• * 


• * 


U^L (k) b 0* 

_ IF (Ov (K))300f JUb f 300 

' 3o6 ns s ,iS ♦ 1 

OLMOX B abs (olm • ox (K)) 

UtL («) » CLOSEN ( - OLMW, ,TNA (NS, Ni|M ♦ I>» DL HQX) 

XTNA (NS. NUH ♦'ir « del (KJ 

30b COnTImUE 

535 COnTImUL 

ZZ * 7 

3)0 IF (N .UT. 1) GO TO b^O 

do 31 k K B 1. N 

315" X <X) s X (K) ♦ ZZ * DEE IK) 

540 CONTI.'UE 

Call talc «x. el. vn) 

A s 0 . 

B B 0. 

C B 0 . ; 

D a 0 . 

IF (M .LT. 1) GO TO b4B 

do 32 k 1 a !• M 

DE~(rr»*0. 

IF (Ni.M .lT. 1) GO TO bbo 

do 32^ X B 1. NUM 
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i 

L 





I»E U)- ♦ XTM (It K) -*-xTNA (Kt-NJM ♦ 1) 

•OE (I)) / Z •• 2 


• E (I) • nOE <I) 


3 iU-Ot 11 ) • 

550 CONTlrUt 

Due (t» « (Efc (1) - t (I) - 

A B A-*'0E "(1) • OOE 

d B B « OE (1) * (I) * 2 

C B C ♦ E (I) *01 (I>‘ 

'325~0 -« 0 ♦ DDE 
545 continue 

A B 3.0 • A 

D B 0 ♦ 0 

ZP B . C / a 
DO 33" 1 B 1* 10 

OZ B -,C ♦ ZP 0 1b~ > I P • - ( A ~ZP^~DTTT~/~T8~»~Z^ 

IZP • r) ) 

Zp B 7P - DZ 


» - o - 


330 

335 


'340 

55a 

345' 

350 

355 


360 

560 


IF (A*.S (OZJ -•' 0 , 00011335 ,— 3 ^ 5.--330 
CONTIttOt 

ZP b rLOSEP <ZM* ZP-, i.'J) 

6 QE » 0 .* 

If (m ,lt. n GO TO 555 

bo 34 " I B !. t M . _ 

iCE * SUE ♦ <E iri-*'ZP • (DE"'tl) ♦ 4 OOC"IPII • • I 

CONTInUt 

VP B sQHT (SUE / Fm) 

BET B (V - VN)- / V • 100 .“ 

IF (8rT)370f 370, J55 

V B Vm ' 

KX ( 8 ) « 1 ~ ‘ 

IF (M ,UT. l» 60 TO 560 
00 36r I B 1* M 
E (1> » EE (1) 

CdNTIr Ut 


' 365* 


370 


3b0 

385 

390 


If (drT - PEH)3f5. 365 » 365 

Z B Zp 
KX (7, » 1 
KX6 3 KAb - 1 
IF (Kvo)410» 410, 200 
KX6 s KXb - AX U) 

IF (Kyd)430, ‘♦30, 440 
375 — KX 6 » KX6 - AX I 2 ) 

IF (Rv6)4l0, 410, JdO 

IF (1«0. * (V -*VP> / V - PE'<)410, 385, 385 
ZZ 3 7 P - Z 
IF (Z7>390« ‘»10* 405 
ZZ 3 . ZZ 

IF (N ,LT. l» 60 TO 565 
DO 39e A s If N 
395 OtL (f ) * “ 

- 565 “ Continue 






IF (NiiM ,uT. 1) 60 TO 570 
bo 40" A 3 I, NOM 

40o XTNA (A, MJM't 1) b"-“XTNA (Kf NUM *1) 
570 CONTIf.Ufc 
405 Z B Z 7 

GO'TO 310 

410 Z B Zo 

AX (Ir) * 1000, * <00 - VJ / VO 
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IF (K» (10) “ KA“(6) 415 ~ ■ 

41S Z « Zp 

420 (F (N .LT. 1) GO 

“bo 42c L ■ It N 
425 Y (L) ■ X (L) 

575 COisiTU-Ut 

U a V 

W a Z 

HETUMm — 

•4J0~ IF (r» .l.m)~GO T0~5B0 
UO 43e I a It N 

435 A (!)_■ X (I> - Z • U^L (I) 

■5bo~C0NTIv0e 

60 TO 410 , ^ 

440 IF tlrO« * (V - VP) / V - PE >430t 44S, 4 45 

"445 2Z ■ 7^ - i 
Z a Zn 
60 TO 310 


-SUbROiiT iNE <»AoSitL- -<Cf-NRO*~^Rr NCC* 

UlMENelON C (15» i&)t L (16#‘ 2) 

bITS ; 2. • * ( - lb) 

'NH » 

NC > i^ CC 

4NITULUE. 

N5F ■ 0 
NhiM ■ NH - I 

NHP » N« ♦ I 

-0 . 1. - 
LSO s 1 

IF (Nn .LT. 1) bC* TO 22S 
"00 10» KR B-i , NR ■■ 

L (KR, U ■ ^R ’ 

L (KR. 2) ■ 0 

"COnTU.UE 

(Nr - l)l05t Ibbt lOb 
ION RHAbt. 

IF (NpM- .LT. l)"0(J-TCrW9 

do ibA KP * 1 * NHM 

KPP B KP ♦ 1 ■ _ 

"PM ■ r; ■ 

HPN * 0 

SEaWCu column KP from UIAGONaL down, for max pivot. 

-If (no .lt.-kp) 

do 11c Kh » KP* nr • 

LKR B L (kR* 1) 

"PT B .Ub (C (LKH»"KP^ 

IF (Pt > PMlilb* iiS* llo 

Pm b r>T 

• MPn b kh 
LMP B L'XH 

CONTlrUfc 

coNTi.ut; 

IF MAv PIVOT IS ZERO* matrix IS SINGULAR* 

IF (MpN)l?0« 2lb* l<iV 

NSF » N5F ♦ I ' ‘ ■ 

IF (MpN - KP»12b* IJO* 125 

NEW HpW NUMoER KP HAS MAX PIvOT. 

LSO » - LSD ■-- - 
L (MP... 1) * L IKP. 1) 

L (KP: 2) b l (MPNi U 

L (KP, 1) B LMP ■ ■ 

HOW Operations io zero column kp below diagonal. 

MKP B L (KP* 1) • 

P B C (MKP* KP) 

() B 0 • P 

IF (Np .LT. KPP) go 10 240 
— do 15J KH B KPP* NR ■ ' 

MKR b L (KR* 1) 

(J B c (MKR* KP) / P 

IF (0) Ub. ISO* I3b 

SUBTR.CT 0 * PIVOT ROW FROM r.Ob KR. 

IF (Nr .LT. KPP) do TO 245 | 

bo I4c LC « KPP* NC 
R B 0 * C (MKP* LC) 

C (MKo* LC) * C (MKH* LC) - .» 



IF (A..S (C LO) - ARS («» • BlTS)lA0t U5. 1*5 

I4U C «Mm.* LC) ■ 0 . 

Us tONllnUt 
?4»3 COMIfUt 
ISU CONllrUL 


<f*u 

23U 

C 

155 

CUriTIfUt 




C0MI»iUt 

LOwFH HlGhT rIANU CoRNtH. 
LNR a L (NH» 1) 



160 

P a c (LNHt 'XH) 

IF (P)loo. 215* l60 

NSF a N5F ♦ 1 



c 

U a D • P * L5U 

IF (No - NO 165* 210* lo5 

bACK SOLUTION pMASt. 

• 


“165“ 

IF (Nr .LT. NRPrGOlO 250 
bo l9r HC ■ NRPt NC 
C ILNr # MC) ■ C (LNH* MC) / ? 

• 


170 

IF (Nr>~- 1) 170* 190*~IJ5 
IF (NoM .lT . U GO TO 265 
bo Ibe LL « 1* 


• • 


KR a ViR “ LL 
HR a 1 (KR* 1) 
KRP a kH « 1 




IF (No .LT. KRP) 60 TO *260 
bo I8r rtS a KRP* NH 
LHS a L (**S* 1) 







R a C (MRt MS) • C (LMS*“MCT 
C (MR, MC) “*C <MR, Mc> - R 
IF (ArS (C (MRt HC)) - ABS (i(), 

* HITS) 175* lao* 180 


■ 17s“ 
ISO 
260 
166 
255 
190 
2S0‘ 
C 

C (MR. MC) * 0^^ 

Contu.ol 

CONiT T. UlE. 

• 

*• 

C (MR, MC) * C (MR» MC) / C MR. KR) 

CoNTI.’ut 

(iONTlnUt 

#« 

■ COnTI* Ut 

ShuFFi t SOLUnOM KORS bACK 

If (N,rM .lt.'i) go to 205 

natural OpOER* 

«• 

bo 20c LL » 1* 
m - LL 

MKH a L (KR* ?) 

196 

IF (M.,H)ig5« 20?* 195 

MKP a L (KR* n 

IF INr- .LT. •'<RP) GO 10 270 


• • 


bo 20;i LC * nrP* NC 
0 a C (MkR* LC) 

C (MKp, LC) » c (MKP* LC) 



2 oo 

270 

205 

' C (MKP* LC) ■ Q 
conti-ol 
CONTImUL 


• • 

- 265 
C 

210 

COnTImUE 

NORMAi AND singular MLIORnS. 
C (1, 1) a 0 

6000 SOLUTION COULD HAVE O»0. 

• 

215 

220 

-GO TO 220 
C (1. 1) a 0. 

Aetur*. 

• 



115 

-^0 
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SOhROmTINP SOKftSS^H* 

OlMENclON SS(20)' 

Mo»n 

2 lF(MO.lS) 2 l,^ 17 f< 

IF 

22 MOs2*|HU/4,)*1 

' 6o To 2^" 

23 M0«2«,MU/rt)*l 

24 KOmN»mU 

JO«l 

25 I»JO 

26 lP«l*nO 

27 YT-SsIp 

SS(P;SS(IP) 

SS(IPi*TT 

1F(I-T)2tt«26f26 

28 J 0 «j 07 I“ 

If ( JO^KO)25*i^Stitf 

9 HETUHi^i 
tUQ 


ORIGINAB PAGE IS 
OF POOR QUAinY 
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-SUBROiiTiNC t»!»»LOl tPHt T, N» HC) 

COHMOm /PaRAMS /IU» AMIN, XM(>X, YMlNt YmAX« NCYC 
“COMMO^i/MEANS/PBAHUbJ t fUARdB) » 1 HEA 0 ( 8 ) 

Ul«ENc;lUN ID < 6 )» XS ( 2 ). YS « 2 )t XB ( 2 M YB < 2 ) 
"OlMENcION NXH * 16 ) t Nxr ( 6 ) ♦ MYT" (677 N»^T“t 6 ) 
OlMENcIUN XN ( 2 U 0 )« YN ( 200 ) 

UIhEN<ION PP ( 2 Y 0 )t i ( 200 ) 

■p I mEN« I ON* I VCRTlSr^ 1 MOR I Z ( 3 ) 


- 0ATA>XS 
UATA fXo 
UATA (NXh 

-Data *iNXT 

1*1. 0. -it 

UATA(tVERT« 1UM 
— UATArfMUBZ*l(W 
£ntRY TlNlT 


165;~^19;> ♦“( YS ■165, =~T5*.“» 

O*. 160.)'~(YB • •«* ISO,) 

s 1. if 0. 0. 0. 1) ... _ 

. e,--m)rO»' 0 . 2) *^ 'TNvr ir- 3 — r.—i 7 o# o* 2 )» -(nht 

0. 2) 


-CALL" 

Call 

Call 

Call 

Call 


P*YH,10Hr-SSUR£*l (*,10H9MB*1) 

rrvciiK i I OKFKA TORE* i"rrro«»Y<^tiO" 


oOROEi' 

-YSf”or- lOf^ — y~ 
Yb* 1. 1') 


rRfPL'T ■ (XS* 

rHiPLT (X«t 

«EAOI 6$ 

rHiPLT t --IS,.- 102,5. -XT. iHEAtrr^ror 
rSipLT ( - IS., sti.. NVt* IVERT. l«) 

rRfPLT <50.. - 1S;7 NHT. IMO«Z. lO) 




0..-1. ir IT ' 
ISO*, i. 1. 1) 


100 

c 

105 

110 

115 

120 


AH ■ T5u, / (XMAX - XMlN) 
dH ■ . am * XMIN 

call rRIPLt (0* H-3tl-Y-7"*3.“^)' 

Do lOn i B 1. Y 
XX * "5 * I 
call rRTPLT (XX. 

Call rRlPLT (xxi 

CO.'^TImUL 

- YTICWARKS* 

NVPATw ■ 1 * 

If (NrYC ,EQ« 0> NvPATm b 2 

GO TO (105. 110) ♦ NVPAfTI 

B AlOGIO (YMAX) 

B AlOGIO (YMIN) 

115 ^ 

B ymAX 
B YHIN 

750, / (XVMAX - XVrtlN) 


XVHIN 

(O.'O. 


XVHAX 
XVMJN 
GO TO 
XVMAX 
XVMIN 
AV B 

dv B . AV • 

Call rMfPLT 

P b"Ym1N - 

60 TO (120. 

OP » o • .1 
IF -{NrYC ,LT. D GO 
bo 13o NC B i. fiCYC 
00 12e 1 B 1. 9 
-ir— rxf iEO.-V7- 

P B P -'dp 

PV B aLOGIO (p> 


JH*0“. 3. 5) 


135). NVPAtM 

TU-235" 


NCYC) ) GO TO 13ir 



i 


!J 



i' 

i 


j 






YY ■ fV • py ♦ ov — - • 

r)^*HLT (U« YYi It It 1) 
CALL rHiMLT OSUtt YVt It It 1) 



12» 

fcONTltUt 

M ■ l)i> 

UP ■ r. • ,l 
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239 

CONTiV'Ut - - - 

CONTU'UL 

Go TO US 


#• 

H5~ 

OP *'T00. 

DO 14^ 1 ■ It 6 
P ■ P - OP 



UO 

YY • dV 8V 

Call rRlPLT (ot YYt It It 1) 

call rrtfPLT (IbO.t YYt 1» It 1) 



c 

US 

' x-scA,T "vjuon 

OX ■ |XMAX • XMlN) • *2 
00 is: I • 1 t 6 



XV « vMlN • -^T »^x 

XX > to • (I - n - 
^NCOUr (5. 15U0> lin'XV 

^500 

150 

FOMHAt (FS.0> 

Call rHIPLT (XXt - 9tt NXHt lilt lO) 

Continue 



C 

ISb 

r~ s^alE Valued 

60 TO (ISSt 165> t NVPAfH 
OY « TSu, / ncyC 




NCP » UCYC * 1 
P B YmIiY * lu. 

IF (NrP .lT. 1) 60 TO 2<*0 


## 


"OO 1 6 - I B • T t~ NCP 

P » P * .1 
YV B ‘LOGIO (P) 



1502 

Yv a >V • YV * bV 
£nC00f (4t 1802* IIIIP 
FOhmAt (F6.1) 



uo 

240 

Call r8IPLT ( - l2.t~YVft nXM, IHt lo> 

CONTUUt 

COnT ImUE 


• • 

169 

RETURk. 

OY B ,-00. 
P B IrOJ. 



UO 17n 1 s li 8 

YV ■ o - (1 - 1) • OY 

Yv » *V • YY ♦ bV 

170 

ENCOUf (St 1900t IIUYY 

call rRfPLT ( - 10. t 'YVt nXM, lilt 10) 

RETURm 



C 

ENTRY TPLOT 
M B N ♦ 1 



175 

60 TO (17St lB0)t MC 

Call fRIPlT (Ot ’ w. jh.ox, i, s) 

So TO US 



1 8 or 

189 

Call rRiRLT^cov Ot SKvoo^r-it si 
If (h ,lt, 1) 60‘To 245 

bo 23r I ■ It M 


• • 


i 

1 


i 

i 




I 







i 




i 


i 


i 


.jMihMtM Ml III IMlIliSiiiiig 


if II .60. 

190 PVAL * (1) 

TVAL .III) 

— “do TO 210 
19b IF (I .to. M»20Ut 20b 

200 i»VAC ;• HP II - 1) 

i 

60 TO 2i0 

20a ^VAL ; PHAR (I - 1) 

TWAL*;'TBAir IT -'TT 

210 Xn (It ■ AH • TWAL * BH 

GO TO I215» 220) I NVPMiH 

—215 — YN II> *‘AV'^"A1.0G1V'^^»l) . ®V 
GO TO 22S 

220 TN lit » AV • PVAL ♦ bV 

— 226 -Call rRiPLT-txN-nr» t.N m* li t» ir 

230 SONTU'Ut 

24b COnTW-UIl _ 

Call~>r iPLT"tot-ir« o» Of— ti : 

Call rRiPLT (xNt yn» h» i, ij 

RETURm 

4nd — 


